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Table for Flow Conversion 
























































Unit m/sec | m/’/day Usec ft*/sec fi'/day | ac-idday | gal/min | gal/day mgd 

1 cubic 

ome 1 964x10°| 10 35.31 |3.051x10°| 70.05 | 1.58x10* | 2.282x10"| 22.824 
second 

1 cubic 

meter) | 1157x105} 1 0.0116 |409x10*| 3531 | 81x10*} 0.1835 | 26417 |2.64x104 
day 

| liter 0.001 86.4 1 0.0353 | 30512 | 0070 | 1585 | 2.28x10*| 228x107 
second 

1cufooY | 60283 | 24466 | 2832 1 |86¢x10'| 1984 | 4488 | 646x10°| 0.646 
second 

1 cu. foo 

‘- 3.28x107| 0.02832 | 3.28x10| 116x105} 1 23x10° |5.19x10°| 748 | 7.48x10° 
1 acre-foo 

“= 0.0143 | 12335 | 14276 | 05042 | 43,560 1 226.28 | 3.259x10°| 0.3258 
I gallo | 63x10°| 5451 | 0.0631 |223x10°| 1925 | 442x10°] 1 1440 «| 1.44x10° 
. ow 4.3x10* | 3.79x 10° | 4.382x 10°) 1.55x10*| 0.11337 | 3.07x 10° | 6.94x 10* 1 10¢ 

1 million 

gallons) | 438x107] 3785 43.82 155 1133710] 3.07 694 10° 1 

day 











Conversion Vaiues for Hydraulic Conductivity 








1 gal/day/ft’ 
1 gal/day/ft? 
1 gal/day/ft? 


1 fUday 
1 fUday 
1 fl/day 


1 cm/sec 
1 cm/sec 
1 cm/sec 
1 m/day 
1 m/day 
1 m/day 


0.0408 m/day 
0.134 fi/day 

4.72 x 10° cm/sec 
0.305 m/day 

7.48 gal/day/ft? 
3.53 x 10* cm/sec 
“64 m/day 

2835 fi/day 
21,200 gal/day/ft’ 
24.5 gal/day/ft’ 
3.28 fi/day 
0.00116 cm/sec 
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Abstract 


This guide analyzes the field methods involved in conducting a geohydrologic analy- 
sis, including pretest water level monitoring, pumping phase, and recovery phase. 
Selected methods of analytical analysis are reviewed with reference to tlh geohy- 
drologic setting, the stress placed on the aquifer by the pumping well, the observa- 
tion of aquifer response, the mathematical solution to the hydraulic head response 
in the aquifer, and the technique for calculating the hydraulic properties of the 
aquifer. Type curves are included for selected aquifer test methods. 
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Introduction 


Purpose and Scope 


The need for a comprehensive geohydrologic analytical guide for Bureau 
of Land Management field offices became apparent as geohydrologists 
and water resource specialists were called upon to interpret and evalu- 
ate data in support of ground water resource projects, with specific 
emphasis on mine dewatering projects. Today’s mining operations take 
place, for the most part, in the form of open pit and underground work- 
ings, all of which require, to some degree, dewatering of geological mate- 


rials for mineral extrav:tion anc for safety. This is particularly important 


because of the increasing emphasis placed on water resources nation- 
wide. The purpose of this guide is to provide methods for geohydrologic 
analysis as applied by geohydrologists and water resource specialists 
working on mine dewatering and other water resource projects. 


The guide presents analyses of a variety of ground water problems en- 
countered in the planning and development of Environmental Assess- 
ments (EAs) and Environmental Impact Statements (EISs) for mine 
include analysis of depletions caused by pumping, estimated seepage, 
analysis of drawdown, and estimates of permeabilities for hydrostrati- 
graphic units. In analyzing these and other ground water problems, 
theoretical assumptions and limitations are outlined and specific meth- 
ods are addressed through the use of tables, figures, and solution equa- 
tions. 


Previous Work 


An extensive Summary of Hydraulic Test Methods is presented in Ap- 
pendix A of this work. Of the original papers describing hydraulic test 
methods, the paper by Theis (1935) is highly recommended reading for 
the interested geohydrologist or water resource specialist. Theis’ (1935) 
paper introduced the most useful method of aquifer flow hydraulics and 
aquifer flow concepts. In addition to these original papers, the 
geohydrologist or water resource specialist will find the Selected Refer- 
ences section useful in providing guidance on selection of aquifer test 
methods, interpretation of aquifer test data, and examples of applica- 
tions of hydraulic test methods. The report by Stallman (1971) is useful 


for the practical application of aquifer test planning and data interpreta- 


tion. The report by Lohman (1972) is an extremely good text on the basic 
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principles of ground water hydraulics and methods with examples of 
their application. Reed (1980) presents the most complete collection of 
tables and types of curves for application of aquifer test methods io 
confined aquifer problems together with discussions of the «nalytical 
solutions and their limitations and applications. Other use '\)! references 
include Ferris and others (1962), Walton (1962), Bentall (1962a), 
Hantush (1964a), Kruseman and DeRidder (1991), Dawson and Istok 
(1991), Fetter (1994), and Vukovic and Soro (1991). Walton (1962) gives 
many examples of aquifer tests including information on the geohydro- 
logic setting, test data, and type curve applications. These references 
from the early 1960s are outstanding treatments of many useful hydrau- 
lic test methods. In addition, several important methods have been 
developed in recent years, such as methods for unconfined aquifers, 
pumping well storage, inertial effects, advances in slug test procedures, 
and solutions to boundary value problems by numerical inversion tech- 
niques (Moench and Ogata, 1984). 





Geohydrologic Analytical Procedures 


Geohydrologic Characteristics Determined From 
Aquifer Tests 


An aquifer test is a controlled in situ experiment made to determine the 
geohydrologic characteristics of water flow and associated rocks. The 
test is made by measuring ground water flow or head that is produced 
by known hydraulic boundary conditions such as pumping wells, re- 
charging wells, variations in head along a connected stream, or changes 
in weight imposed on the land surface. 


The geohydrologic characteristics that can be determined from an aqui- 
fer test depend on the onsite test conditicns and installations. The most 
common geohydrologic parameters determined are the coefficients of 
transmissivity, T, and storage, S, or storativity. Transmissivity is a 
measure of the ease in which the full thickness of the aquifer transmits 
water; the hydraulic conductivity, K, is a measure of the ease with which 
a unit thickness of the aquifer transmits water.* Therefore, 


where b is the thickness of the aquifer. The evaluation of rates of 
ground water flow in an aquifer requires knowledge of hydraulic conduc- 
tivity and effective porosity. Effective porosity, or d ainable porosity, is a 
measure of the interconnected void space of a medium. The storage 
coefficient of an unconfined aquifer is approximately equal to the effec- 
tive porosity, n. The storage coefficient of a confined aquifer is typically 
much smaller than that of an unconfined aquifer. Whereas water yielded 
to a weil ffom an unconfined aquifer is derived principally from drainage 
of water from voids, water yielded to a well from a confined aquifer is 
derived principally by compression of the aquifer and expansion of the 
water. Values of effective porosity of granular materials usually range 
from 0.1 to 0.4; storage coefficients of confined aquifers usually range 
from 10* to 10°. 


The relation of flow in an aquifer to the hydraulic conductivity and the 
hydraulic gradient are expressed by a general form of Darcy’s law: 


ve-xkKG@.2 (2) 





* For the definition of the symbols used in a specific equation, the reader is referred to 
Appendix B, List of Nomenclature and Symbols. 
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where v is the flux specific discharge, also called Darcy’s velocity, dh/dl is 
the hydraulic gradient, Q is the discharge, and A is the cross-sectional 
area. From these parameters of the aquifer, the rate of advective trans- 
port of a solute can be calculated by the following relation: 


v= (X) (& 
v= (=) (5) (3) 


Thus, aquifer tests do not provide a direct analysis of the parameters K 
and n. But, K can be determined from an aquifer test where the satu- 
rated thickness is known. The effective porosity can be estimated as the 
storage coefficient from tests of an unconfined aquifer. The determina- 
tion of storage coefficient from an aquifer test requires analysis of the 
drawdown response in observation wells rather than in the pumping 
well. Drawdown response solely in the pumping well can be used to 
calculate transmissivity, but is not reliable for determination of the 
storage coefficient because the effective radius of the pumping well is not 
known (after Bedinger and others, 1988). 


Application of Hydraulic Tests for Geohydrologic Systems 


Aquifer tests were originally applied to wells completed in aquifers that 
were used for water supply. The first tests were designed simply to 
define the gross hydraulic properties of the water-yielding material. The 
earliest application of aquifer tests was in the design of well fields and in 
the prediction of the performance of an aquifer as a source of water 
supply. Aquifer test methodology has increased tremendously in sophisti- 
cation as a result of more complex techniques applied to analyzing 
simple aquifer and boundary conditions. The type of aquifer test meth- 
ods available today can provide more detailed information on the confin- 
ing beds as well as flow system characteristics. Aquifer test methods 

can provide much more of the detail needed for characterization and 
analysis of hydrologic systems. 


Definition of hydraulic properties is an essential element in character- 
ization of geohydrologic systems and in design of ground water field 
programs for mine dewatering and water resource studies. Aquifer test 
methods are specifically designed to provide analysis of hydraulic prop- 
erties under a certain set of geohydrologic conditions. Therefore, aquifer 
tests can be designed and performed to provide the type of information 
on the flow system that is best suited for the geohydrologic setting and 
the application for which the data are needed. As discussed by Bedinger 
and others (1988), aquifer tests can be chosen to provide information on 
the gross hydraulic properties of a large volume of an aquifer, the 
hydraulic conductivity of a relatively thin bed in the flow system, the 








Geohydrologic Analytical Procedures 


relative horizontal to vertical hydraulic conductivity of an aquifer, areal 
anisotropy of the aquifer, or the leakage from a confining bed. 


The aquifer test method chosen must provide the type of information 
required for a given application. For example, pit and underground mine 
operational monitoring and mitigation programs might be enhanced by 
information that includes horizontal and vertical hydraulic conductivity, 
estimation of the rate and direction of ground water flow, spatial distribu- 
tion, and the hydraulic characteristics of a specific hydrostratigraphic 
unit. In addition, the design of a plan for ground water reinjection or 
infiltration for mine operational water management might require one or 
more long term aquifer tests with many observation wells. 


Aquifer tests require information on the geohydrology of the area and a 
network of one or more wells that are constructed and instrumented to 
provide the data necessary for analysis by the aquifer test method cho- 
sen. Unless the test area has been defined by investigations such as 
borings, geophysical logging, coring, surface water surveys, water level 
measurements, or other means, the most appropriate aquifer test method 
may not be chosen. Aquifer tests designed for analysis of specific hydrau- 
lic properties generally have specific requirements for layout and con- 
struction of the pumping and observation wells. 


Hydraulic Test Planning, Design, and Implementation 


An outline of the steps involved in the planning, design, and implementa- 
tion of an aquifer test is given in the following sections. 


Evaluation of the Geohydrologic System 


Through an evaluation of the geohydrology of a water resources project, a 
conceptual model of the flow system is made. The evaluation needs to 
provide a concept of the nature of the aquifer’s transmissivity, homogene- 
ity, and isotropy, and whether the aquifer is confined or unconfined and, 
if confined, whether the aquifer is overlain or underlain by leaky or 
nonleaky confining beds. Emphasis is placed on the need for knowledge 
of the geohydrologic setting because the response of an aquifer system to 
stress is not unique to the geohydrology. Misunderstanding the geohydro- 
logic setting could lead to selection of an inappropriate aquifer test 
method and incorrect analysis of hydraulic properties. Therefore, an 
accurate estimate of the flow system characteristics needs to be made, 
and an estimate of the hydraulic properties of the aquifer needs to be 
known in order to plan and implement the most representative aquifer 
test. 
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Survey of Selected Aquifer Test Methods 


The available literature on aquifer test methods is extensive and each 
method is specific with respect to geohydrologic conditions and well 
control. Furthermore, each method is usually limited to a relatively 
simple set of aquifer characteristics and boundary conditions as opposed 
to the complexity of the actual area being studied. Selection of the aqui- 
fer test method as discussed by Bedinger and others (1988) is made on 
the basis of the geohydrology of the test site and the field test conditions. 
The geohydroiogy of the test location with regard to nonleaky confined 
aquifer, leaky confined aquifer, unconfined aquifer, and other natural 
conditions of the area determine the applicable set of aquifer test meth- 
ods. The field test conditions (with regard to number and location of 
observation wells, if any), instrumentation for measuring water ievels, 
screened interval, and capacity of the pump on the pumping well, deter- 
mine which aquifer test methods can be applied to the data. These and 
other factors determine the physical constraints on stressing the aquifer 
and on determining the aquifer response, and may further limit the 
aquifer test methods applicable for analysis. 


An overview of some of the more commonly used aquifer test methods 
and their applicability to geohydrologic conditions and field test condi- 
tions is provided in Table 1. Each test is discussed in the Hydraulic Test 
Methods for Aquifers section with information on the applicability of the 
methods to specific test site conditions. General guidelines for the num- 
ber of observation wells and the distance of observation wells from the 
pumping well for the aquifer test methods are given in the next section. 


Well Siting and Screened Intervals 


A single well test uses the same well as the pumping well and the obser- 
vation well. Many other aquifer test methods can be applied to the data 
from the pumping well. The applicability of the common aquifer test 
methods are outlined in Table 1. Determination of the transmissivity is 
considered representative by single-well test data, but determination of 
the storage coefficient is considered unreliable because of the problem of 
estimating the effective radius of the pumping well. Slug tests are 
commonly conducted in wells screened through only part of a saturated, 
permeable section. Tests in such wells measure the properties of only a 
small part of the water-yielding section. These measurements can be a 
benefit when information on variations in the hydraulic conductivity at 
many points is desired. 


The distance from the observation well to the pumping well, r, will 
usually be discussed with reference to a distance, 








Table 1. Aquifer test methods 
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<1.5b,| 4% 
De 1.sb | (4) 


where b is the aquifer thickness, K, is the vertical hydraulic conductivity 
of the aquifer, and K, is the horizontal hydraulic conductivity of the 
aquifer. Obviously, K, and K, are not known prior to a test, but they can 
be determined by a few tests. A general rule of thumb used by many 
geohydrologists where K,/K, is not known is to estimate Do as 2 or 3 
times the thickness of the aquifer. For a fully penetrating pumping well, 
observation wells can be fully or partially penetrating and either within 
or without a distance Do = 1.5b (K,/K,) from the pumping well. If the 
pumping well partially penetrates the aquifer, the observation wells can 
be either fully penetrating within a distance of Dp from the pumping well 
or they can be partially or fully penetrating outside a distance of Do from 
the pumping well. No observation wells are used in slug tests; the slug 
well should be fuily penetrating, but usually is partially penetrating. In 
applying the radial-vertical methods of Hantush (1966a and b) and 
Weeks (1969) to determine horizontal and vertical permeability, the 
pumping well needs to be partially penetrating; the observation wells 
need to be piezometers that are either open at a point or screened for only 
a short vertical distance. The piezometers need to be within the distance 
Do of the pumping well. In applying the general method of Neuman 
(1975) for unconfined aquifers, the fully penetrating pumping and obser- 
vation wells must meet the requirements of distance Do from the pump- 
ing well. The family of type curves for this method is presented in the 
Solution of Boulton and Neuman section of this guide. The method of 
Neuman (1975) requires fully or partially penetrating wells with greater 
or lesser distances of Do from pumping wells to observation wells. 


The nonequilibrium method of Theis (1935) is applicable to unconfined 
aquifers where the pumping well is fully penetrating, the observation 
well is fully or partially penetrating, and the observation well is greater 
than b/r(K,/Ky) from the pumping well for times greater than 10Syr’/T 
(Neuman, 1975, p. 337), and where S, is the specific yield. There are two 
zones where this method can be applied using fully penetrating observa- 
tion wells. The first zone is far from the pumping well at later times 
where r > b/r(K;/K,) and t > Syr’/T (Neuman, 1975, p. 338). The second 
zone is near the pumping well at early times where r < 0.03 b/r(K,/K,) 
and t < Sr’/T (Neuman, 1975), where S is the storage coefficient. 





Establishing Baseline Water Level Fluctuations 
Water levels continually fluctuate in response to local or regional stresses 


that are imposed on the flow system, such as recharge, discharge, 
changes in hydraulic head at the boundaries of the flow system, 
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barometric changes, and weight on the land surface. The effects of these 
background water level fluctuations in the locality of an aquifer test 
ideally are smail; but even so, they cannot be discounted. 


Measurements made before and after the aquifer test to detect regional 
water level trends are required to interpret the background water level 
trend and to more accurately identify drawdown. The water level trend 
before and after an aquifer test at a theoretical observation well is 
shown in Figure 1. The effect of drawdown imposed by the pumping 
weil is superposed on the background water level fluctuations. The 
drawdown is the distance between the background water level trend and 
the water level in the cbservation well. From the theory of ground water 
hydraulics, it is noted that the recovery period is longer than the pump- 
ing period. The recovery of water level after stoppage of pumping is 
measured from the interpreted drawdown curve. 


An inverse relation between barometric pressure and change in water 
level in confined aquifers is commonly identified. Water levels in un- 
confined aquifers are unaffected by changes in barometric pressure. 
Water levels in confined aquifers need to be corrected for barometric 
changes during an aquifer test according to the barometric efficiency of 
the aquifer. 


Step Drawdown Test 


The step drawdown test is usually conducted to provide a basis for se- 
lecting the discharge rate for a long-term aquifer test. A step drawdown 
test is a preliminary aquifer test that uses incremental increases in the 
pumping rate starting from an initial slow pumping rate to successively 
faster pumping rates. The test usually is conducted in 1 day. Pumping 
times need io be the same for each rate; either the water level may be 
allowed to recover between pumping periods or the pumping rate may be 
increased without a recovery period. The duration of each step needs to 
be long enough (usually 1 to 2 hours is adequate) for the rate of draw- 
down to become virtually stable (Figure 2). 


In a pumping well, the major part of the drawdown occurs in the forma- 
tion where the energy provided in overcoming the frictional resistance of 
the formation against the slowly moving water is directly proportional to 
the rate of motion. Another important part of the loss is a function of 
the proportionality of the velocity approaching the square of the velocity. 
A relation between these two components of drawdown is expressed by 
Jacob (1947): 
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Figure 1. 
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Hydrograph of hypothetical observation well showing background 
water levels, and drawdown and recovery of water level. Drawdown 
begins at t=0 and ends at t=1. Recovery is not complete at t=2. 
Residual drawdown at t=2 e..uals the drawdown that would have 
occurred from t=1 to t=2 had pumping continued at a constant rate 
(After Vukovic and Soro, 1991). 




































































Figure 2. Hydrograph of a step-drawdown test: t is equally spaced time(s), s is 
drawdown, and Q is the discharge rate(s) (After Vukovic and Soro, 
1991). 
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8, = BQ + CO (5) 


where 8y is the drawdown in the pumping well, B is the coefficient of 
head loss linearly related to the flow, and C is the coefficient of head ioss 
due to turbulent flow in the well, aquifer, and across the well screen 
(Cooley and Cunningham, 1979). Components of drawdown are shown 
in Figure 3. Rorabaugh (1953) presents a more general form of equa- 
tion, substituting n for the exponent 2. 


Equation 5 expresses the well loss component of drawdown in proportion 
to the square of the discharge, Q. Bierschenk (1964) presents a graphi- 
cal method for determining the constants B and C in equation 5. 


Developing an Aquifer Test Plan 


The following guidelines for specifications and tolerances of measure- 
ments for the aquifer test are primarily from a report by Stallman 
(1971). For additional detail and discussion of these items, the reader is 
referred to Stallman (1971) and Driscoll (1986). These items may be 
used as a checklist that includes tasks that need to be done before, dur- 
ing, and after the test. 


1. Pumping well needs to: 
a. Be equipped with reliable power, pump, and discharge control 
equipment to maintain the discharge rate during the aquifer 
test. 


b. Be equipped to carry discharge water away from pumping and 
observation wells. 


c. Be equipped to measure discharge at specific times during the 
aquifer test. 


d. Be equipped to measure the water level before, during, and 
after the aquifer test. 


e. Have a known diameter, depth, and screened interval(s). 


f. Have a screened interval(s) compatible with the aquifer test 
method. 


g. Be used for a step drawdown test to determine the discharge 
rate for the aquifer test. 


— 
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2. Observation well needs te: 

a. Be used for water level measuremenis during the step draw- 
down test to assure hydraulic connection with the aquifer, 
determine accuracy of water level measurements, and deter- 
mine response to discharge from the pumping well. 

b. Bea known radial distance from the pumping well. 


c. Be used to measure baseline water levels to determine the 
trend of these levels before the aquifer test begins. 


d. Have o known diameter, depth, and screened interval(s). 

e. Have a screened interval(s) and a distance from the pumping 
well compatible with aquifer test methods to be used in the 
analysis. 

3. Aquifer test method(s) need to be: 

a. Selected for analysis based on geohydrologic condition and 
test area installations, especially pumping and observation 
wells and theorized response of flow system. 

b. Known so that applicable type curves, graph paper, and mate- 
rials for onsite analyses of data can be assembled. 


4. Records of and the tolerances in measurements for the following are 
needed for analysis: 


a. Pumping well discharge (+10 percent). 


b. Depth to water in pumping and observation wells below mea- 
suring point (+0.01 ft). 


c. Distance from pumping well to each observation well (+0.5 
percent). 


d. Synchronous time (+1 percent of time since pumping initiated). 
e. Description of measuring points. 
f. Elevation of measuring points (+0.01 ft). 
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Vertical distance between measuring point and land surface 
(+0.1 ft). 


Total depth of all wells (+1 percent). 


Depth and length of screened interval(s) of all wells 
(+1 percent). 


Diameter, casing type, screen type, and method of construc- 
tion of all wells (nominal). 


Location of all wells in plan either relative to land survey net 
or by latitude and longitude (accuracy dependent on indi- 
vidual need). 


5. Measurements of water level need to: 


Be made periodically in all wells 24 to 72 hours before the 
step drawdown test, continuing through recovery (Establish- 
ing Baseline Water Level Fluctuations section and Figure 1). 


Be made continually in all observation wells during the aquifer test. 


Be recorded with a logarithmically decreasing frequency 
during the aquifer test. For example, with discharge com- 
mencing at time zero, measure at 1, 1.2, 1.5, 2, 2.5, 3, 4, 5, 6, 
7, and 8 minutes and at succeeding time multiples. 


Be made continually in all wells after stoppage of pumping to 
determine recovery for a period equal or longer in duration 
than the period of pumping. 


Be made periodically in all wells after complete recovery to 
determine baseline water levels. 


6. Measurements of barometric pressure need to: 


a. Be made continually during tests of confined aquifers, which 


b. 


are affected by barometric changes in water level. Measure 


barometric pressure during pretest through post-test water 
level measurement periods. 


Be recorded to calculate barometric efficiency of aquifer. 


. 
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Analysis of data as it is collected during the drawdown and recovery phase 
is helpful in assessing the progress of the aquifer test and in determining 
the time period necessary for the drawdown and recovery phases. 


The drawdown phase of an aquifer test provides the primary data for 
analysis of aquifer characteristics. Activities that need to be performed 
during the drawdown phase are: 


1. Plot measured discharge and measured depth to water for the 
pumping well and each observation well. 


2. Correct baseline water level fluctuations and drawdown 
water levels for fluctuations in barometric pressure, as 
applicable. 


3. Interpret baseline water level fluctuation from plot of cor- 
rected water level and calculate drawdown (Figure 1). 


4. Plot data for analysis according to the aquifer test method or 
methods selected for analysis. 


5. Evaluate progress of drawdown phase on the basis of analysis 
of the hydraulic properties by the aquifer test method or 
methods selected. This is done by rating the fitting of the data 
to type curves or rating the time at which the data plot is a 
straight line if using the modified nonequilibrium method. 


6. Terminate drawdown phase when analyses indicate that data 
are adequate for calculating hydraulic properties by the aqui- 
fer test method or methods selected. 


The recovery phase provides a data set for several aquifer test methods 
that can be used to verify the drawdown phase calculations. Recovery 
data analyses are considered by some geohydrologists to provide more 
accurate calculations of hydraulic properties. Minor variations in dis- 
charge that may have occurred during the drawdown phase are not 
apparent during the recovery phase. Recovery measurements in the 
pumping well may provide more accurate estimates of hydraulic conduc- 
tivity because well loss is smaller during the later recovery phase. The 
recovery phase provides a transition to the baseline water levels after 
recovery and a basis for re-evaluating the drawdown and recovery water 
levels. Activities that need to be performed during the recovery phase 
are: 
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1. Continue to plot measured depth to water for the pumping 
well and each observation well. 


2. Correct recovery water levels for fluctuations in barometric 
pressure, as applicable. 


3. Interpret baseline water level, interpret plot of drawdown 
from discharge phase, and calculate recovery of water level 


(Figure 1). 


4. Plot recovery data for analysis according to the aquifer test 
method selected for analysis and calculate hydraulic 


properties. 


5. Continue recovery measurements to document post-recovery 
baseline water level. 


After the aquifer test is completed, all data need to be reconsidered and 
revised analyses made as needed. The drawdowns may require revision 
based on final predictions of baseline water levels before and after the 
test. Corrections may be necessary in type curves or drawdowns for 
changes in discharge rate. It may become apparent that aquifer bound- 
aries are reflected in the data and that the effects of such boundaries 
need to be assessed. 


Type Curve Utilization 


The solution to several of the principal aquifer test methods depends on 
the application of type curves to plots of the aquifer test data. The use of 
type curves is required by the existence of integral expressions in the 
analytical solutions that cannot be integrated directly. The application 
of type curves to solve for the aquifer properties follows a similar proce- 
dure in each instance. The type curve solution discussed in this guide is 
for the solution to the Theis (1935) nonequilibrium drawdown method. 
Application of type curve solutions to other methods, for example, the 
methods of Hantush and Jacob (1955) and Hantush (1960) for leaky 
confined aquifers, and the methods of Boulton (1963) and Neuman 


(1972) for unconfined aquifers, follow the same general procedures. 


joo, 
aj 








BLANK PAGE 








Hydraulic Test Methods for Aquifers 


The aquifer test methods discussed in this section are for the simplest 
geohydrologic site conditions. Discharge is assumed to be constant; the 
pumping well is assumed to be a line source. Therefore, wellbore storage 
is ignored, and the aquifer is assumed to be homogeneous, isotropic, and 
areally extensive. Solutions to these principal methods are straightfor- 
ward and type curves are widely available. 


Nonleaky Confined Conditions 


Solutions to flow conditions induced by discharge from a well in a 
nonleaky confined or artesian aquifer are considered first. Though based 
on simple boundary conditions, the solutions to the methods discussed 
here are useful when applied to appropriate geohydrologic conditions. 
The methods may also be applied to obtain a preliminary estimate of 
hydraulic properties as discussed by Bedinger and others (1988) for a 
test in which geohydrologic conditions are not well known, or for a quali- 
tative examination of aquifer test data to aid in selecting an appropriate 
method or model. 


Theis Nonequilibrium 


The solution of Theis (1935) for the change in distribution of head near a 
well being pumped revolutionized aquifer test methodology and the 
study of aquifer hydraulics. Although about 50 years old, Theis’ method 
is still the most widely referenced and applied aquifer test method. The 
Theis solution is the basis and limiting case for solutions to the head 
distribution in many geohydrologic situations. 


Assumptions: 
1. The pumping well discharges at a constant rate, Q. 


2. The pumping well is of infinitesimal diameter and fully pen- 
etrates the aquifer. 


3. The nonleaky confined aquifer is homogeneous, isotropic, and 
areally extensive. 


4. The discharge from the pumping well is derived exclusively 
from storage in the aquifer. 
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Implicit in these assumptions are the conditions of radial flow; that is, 
there are no vertical components of flow and no dewatering of the aqui- 
fer. The geometry of the assumed aquifer and well conditions are shown 
in Figure 4. The Theis (1935) nonequilibrium solution is: 





ul “s e~” 
o- 2. fi Sw 6) 
and 
a 
us= az (7) 
where 
[ <~ dy = Wu) = -0.577216 -log, u + u 
‘ u? ‘ u? < u* ‘ (8) 
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The integral expression in equations 6 and 8 cannot be evaluated ana- 
lytically, but Theis (Wenzel, 1942) devised a graphical procedure to solve 
for the two unknown parameters, transmissivity, T, and storage coeffi- 
cient, S, where 


ss (<2) w(u) (9) 
and 
> 


The graphical procedure is based on the functional relations between 
W(u) and s, and between u and t, or t/r’. 


Steps to perform the Theis procedure are: 


1. A type curve illustrating the values of W(u) versus values of 1/u is 
plotted on logarithmic scale graph paper (Figure 5). This plot is 
referred to as the type curve plot. Values of W(u) for values of 1/u 
from 10" to 9x10" are tabulated by Reed (1980). Values of W(u) 
for values of 1/u from 10°” to 9.9 are tabulated in Ferris and 
others (1962), and in Lohman (1972). 


2. On logarithmic tracing paper of the same scale and size as the 
W(u) versus 1/u curve, values of drawdown, s, are plotted on the 
vertical coordinate versus either time, t, on the horizontal 
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Figure 4. Sections through a pumping wel! in a nonleaky confined aquifer 


(After Fetter, 1988). 
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Figure 5. Theis nonequilibrium type curve of dimensionless drawdown, W(u), as a function of dimensionless time, 
(1/u), for constant discharge from a nonleaky confined aquifer (From Reed, 1980). 
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coordinate if an observation well is used, or versus t/r’ on the 
horizontal coordinate if more than one observation well is used. 
This plot is referred to as the data plot (Figure 6). Alternatively, 
the type curve can be plotted as W(u) versus u and the data plot 
ted as drawdown, s, versus 1/t or r’/t. 


3. The data plot is overlain on the type curve plot and, while the 
coordinate axes of the two plots are held parallel, the data plot is 
shifted to a position that represents the best fit of the aquifer test 
data to the type curve (Figure 6). 


4. An arbitrary point, referred to as the match point (Figure 6), is 


selected anywhere on the overlapping part of the plots and the 
W(u), 1/u, s, and t coordinates of this point are recorded. 


5. Using the coordinates of the point, the transmissivity and storage 
coefficient are determined from the following equations: 


e u 
T oe — ay 


4Tuct 
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Application of the curve-matching procedure to aquifer test data is 
discussed by Lohman (1972). 


Modified Nonequilibrium 
Assumptions: 


The straight line method, also called the modified nonequilibrium 
method, is a solution when u is small and the Theis solution can be 
approximated by the first two terms on the right side of equation 8. 


Solution: 


Cooper and Jacob (1946) and Jacob (1950) recognized that in the series 
of equation 8, the sum of the terms beyond log, u is not significant when 
u = r°S/4Tt becomes small, < about 0.01. The value of u decreases with 
increasing time, t, and decreases as the radial distance, r, decreases. 
Therefore, for large values of t and reasonably small values of r, the 
terms beyond log, u in equation 8 may be neglected. The Theis equation 
can then be written as: 


. ~log, 228 
8 i& [-0.577216 -log, =) (18) 


and 
(12) 
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Figure 6. Relation of W(u), 1/u type curve and s, t data plot (After Stallman, 1971). 
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from which Lohman (1972) derives the following equations: 


— 
T Wake] Blogz s/Biog,,t (14) 


which applies at constant radius and 


T.< rT VIL SCO (15) 
2xAs/Alog,.r 


which applies at constant time. Equation 15 is the same as the Theim 
(1906) equation. 


a 


Equation 14 can be used to determine transmissivity, T, by plotting 
drawdown, s, at a specified distance on the arithmetic scale and time, t, 
on the arithmetic scale versus the distance of the observation wells from 
the pumping well on the logarithmic scale. By choosing the drawdown, 
Ast or 8, to be that which occurs over a log cycle, 


Alog,,t = log,, (@ “1 (16) 
a 
Alog,,r = log,, (=) *1 (17) 
equation 14 then becomes 
» on 
f° a (18) 
and equation 15 then becomes 
oo ote 
Te ~(5ie,. (19) 


The coefficient of storage can be determined from these semilog 
plots of drawdown by a method proposed by Jacob (1950) where 


- 2.30 2.25Tt 
s eal ee (20) 


taking s = 0 at the zero drawdown intercept of the straight line 
semilog plot of time or distance versus drawdown 


(21) 


where either r or t is the value at the zero drawdown intercept. 


and , 
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Theis Recovery 


A useful corollary to the Theis nonequilibrium method was devised by 
Theis (1935) for the analysis of the recovery of the water level in a con- 
trol well. The water level in a pumping well that is shut down after 
being pumped for a known time will recover at a rate that is the inverse 
of the rate of drawdown. The residual drawdown at any instant will be 
the same as if the discharge of the well had been continued and a re- 
charge well of the same flow had been introduced at the same point at 
the instant discharge stopped. The residual drawdown at any time 
during the recovery period is the difference between the measured water 
level and the pumping water level interpreted from the measured trend 
prior to the stoppage of pumping as discussed by Bedinger and others 
(1988). These relationships are shown in Figure 7. The residual draw- 
down, s', at any instant can be expressed as: 


le - - eo” 
come ST w- i, > wn (22) 
The Theis recovery method is applied in much the same way as the 
drawdown method. From the time since pumping ceased t’ becomes 
large, the semilog plot of residual drawdown s’, and the ratio of time 
since pumping started to time since pumping ceased t/t’, is plotted on the 
logarithmic scale. Transmissivity can be calculated from the following 
equation: 


re 2:32 1093, (5) 


4xAs’ (23) 


Slug Test 


The method for estimating transmissivity by injecting a given quantity 
or slug of water into a well was originally described by Ferris and 
Knowles (1954). Included in this category of tests are methods for deter- 
mining transmissivity in a slug well by determining the response to an 
instantaneous change in head in the well. Head change may be induced 
by injection of water, bailing of water from the well, rapid removal of a 
solid cylinder from beneath the water level in a well to which the water 
level was in equilibrium, or application of pressure to the volume of 
water stored in a shut in well. 
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Figure 7. Graph showing drawdown, recovery, and residual drawdown (After 
Dawson and Istok, 1991). 
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Figure 8. Section through a pumping well in which a slug of water is suddenly 
injected (After Reed, 1980). 
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1. A volume of water is injected into or is discharged from the 
slug well instantaneously at t=0. 


2. The slug well is of finite diameter and fully penetrates the 
aquifer. 


3. Flow is radial in the areally extensive, homogeneous, and 
isotropic, nonleaky confined aquifer. 


The geometry of the slug well and aquifer is shown in Figure 8. 
Solution: 

The solution presented by Cooper and others (1967) for wells that are 
not affected by inertially induced, oscillatory water level fluctuations, is 
for a slug well of finite diameter; application of the solution is by match- 
ing of aquifer test data to type curves. The solution and its application 


have been elaborated on by Bredehoeft and Papadopulos (1980), and 
Neuzil (1982). The solution of Cooper and others (1967) is 


« (2%) f° _pu’ ur 
h= (—) I, ({ Exp (5 ) (Jt) 


(uy, (u) - 26; (u)] - Yo (>) (24) 
(uJ, (u) - 2aJ,(u))}/A(u)})} du 


where 





(25) 
r? (26) 


A(u) = [ud,(u) - 2a, (u) }? 
+ (uy, (u) - 2ay, (u) )” (27) 


The head, H, inside the slug well, obtained by substituting r=rw in 
equation 24 is 


2 F(B6,a@) (28) 
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Figure 9. Graph showing ten selected type curves of F (f, a) as a function of B (From Reed, 1980). 
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P(B,0) = (28) f° (axp(-BE) suaiw)) du iggy 
The curves generated from equation 28 are plotted in Figure 9. 


bactitestinn. 


The water level data in the slug well, expressed as a fraction of Hg (that 
is, H/Ho) are plotted versus time, t, on semilogarithmic graph paper of 
the same scale as that of the type curve plot. This data plot is overlain 
on Figure 9 and, while keeping the baselines the same, the data plot is 
shifted horizontally until a match or interpolated fit of the aquifer test 
data to a type curve is made. A match point for B, t, and a is picked on 
the overlapping part of the plots and the coordinates of this point are 
recorded. The transmissivity is calculated from 





(30) 


eer 1 (31) 


As pointed out by Cooper and others (1967), the determination of S by 
this method has questionable reliability because of the similar shape of 
the curve, whereas the determination of T is not as sensitive to choosing 
the correct curve. Figure 9 is plotted from data from two sources (Coo- 
per and others, 1967; and Papadopulos and others, 1973). Tables of the 
F(6,a) are given in Cooper and others (1967) for valves of 8 from 10~ to 
2.15 x 10° and for values of a from 10* to 10°” in order to apply the 
method to formations having a very small storage coefficient. 


Although the method applies to radial flow in a nonleaky confined aqui- 
fer, the method has been applied to partially penetrating \ells where the 
screened interval is much larger than the well radius. In a stratified 
aquifer where the vertical permeability is much smaller than the hori- 
zontal permeability, the flow for a test of short duration can be assumed 
to be virtually radial. The transmissivity thus derived would apply to 
the part of the aquifer in which the well is screened or open. 
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Bredehoeft and Papadopulos (1980) adapted the method for application to 
formations of very low permeability and extended the range of F(8,a). 
Bredehoeft and Papadopulos (1980) described a technique of pressurizing 
a shut in well in low permeability rocks that decreases the response time 
by orders of magnitude. Neuzil (1982) determined that the slug test 
method for very low permeability formations by Bredehoeft and 
Papadopulos (1980) does not assure the condition of approximate equilib- 
rium necessary at the start of the slug test; Neuzil (1982) also determined 
that the compliance of the shut in well and associated piping, which 
determines the response time, can be substantially larger than the com- 
pressibility of water alone. Neuzil (1982) presented a modified procedure 
and testing arrangement for slug tests in low permeability formations. 
The region adjacent to the wellbore may be altered by the addition of 
drilling mud, precipitation of scale, well stimulation, or gravel pack. 
Moench and Hsieh (1985) examined the altered region or skin adjacent to 
the wellbore and concluded that pressure tests may be markedly affected 
by the altered skin. Moench and Hsieh (1985) determined that standard 
methods of analysis are adequate for open well slug tests. 


Airlift Test 


The data analysis procedure is outlined by Kruseman and Ridder (1991). 
This method is very similar to the Cooper and Jacob method (1946), but 
was developed by Aron and Scott (1965) for a well in a confined aquifer, 
with the exception that it is assumed the discharge rate decreases with 
time with the sharpest decrease occurring soon after the start of pumping. 


The procedure involves injecting pressurized air down the well to lift 
water to the surface as shown in Figure 10, while recording drawdown 
and discharge over time. These data are plotted on semi-log paper with 
drawdown divided by production rate (s/Q) plotted on the vertical scale 
and log time plotted on the horizontal scale. This usually results in a 
straight line; the slope of the line is then used to calculate the transmis- 
sivity as shown in Figure 11. A more detailed description of the equip- 
ment setup and layout is presented by Driscoll( 1986). 


The airlift test procedure has been applied on numerous projects as noted 
by Doubek and Beale (1992). Some of the advantages of this method are 
that the test can be conducted with standard exploration drilling equip- 
ment, water level measurements and transmissivities can be obtained, 
and the test is less costly than some other methods. Some disadvantages 
of this method are that the test cannot be used to obtain storativities, 
stresses only affect zones close to the pumping well, and the analytical 
method must meet discharge rate and other constraints as noted by 


Cooper and Jacob (1946) for applicability. 
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Figure 10. Schematic representation of airlift pumping test (Baski, 1978). 














SLANK 















































































































































{ 0 
a 1st Test 
1 ‘, 
v4 
el 
c Bas 
ag 
© 
0 3 2 
3 
1 2 3 48678910 
q Time (minutes) 7 





Figure 11. Airlift pump test drawdown plots (Baski, 1978) 





Hydraulic Test Methods for Aquifers 
Leaky Confined Conditions 


Confining beds above or below the aquifer commonly provide water to the 


aquifer by leakage when the aquifer is pumped. Methods that account for 
leakage will be discussed next. 


Leaky Confining Bed Without Storage 
Assumptions: 
1. Pumping well discharge is at a constant rate, Q. 


2. Pumping well is of infinitesimal diameter and fully penetrates 
the confined aquifer. 


3. Confined aquifer is overlain or underlain everywhere by a 
leaky confining bed having uniform hydraulic conductivity, K’ 
and thickness, b’. 


4. Leaky confining bed is overlain or underlain everywhere by 
an infinite source bed with a constant hydraulic head. 


5. Hydraulic gradient across the leaky confining bed changes 
instantaneously with a change in head in the confined aquifer 
(no release of water from storage in the leaky confining bed). 


6. Flow in the confined aquifer is two climensional and radial in 
the horizontal plane; flow in the ieab:y confining bed is vertical. 


The nonequilibrium technique of Hantush and Jacob (1955), though a 
simplification of a leaky flow system, is widely applied as discussed by 
Bedinger and others (1988). The method assumes an unlimited supply of 
water from the overlying or underlying beds, but no release of water from 


storage in the confining beds. 


The geometry of the assumed well and aquifer system is shown is 
Figure 12. 


The assumption of no release of water from storage in the leaky confining 
bed may usually be met at early times before water is yielded from the 
confining bed and at late times when the system is near steady state. The 
assumption may also estimate conditions for thin confining beds. 
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Figure 12. Section through a pumping well in a leaky confined aquifer without 
storage of water in the leaky confining bed (After Reed, 1980). 





Hydraulic Test Methods for Aquifers 
Solution: 


The solution for the conditions stated as given by Hantush and Jacob 
(1955) are: 





g= 2 f° o-*-(r#/4B*z) a, 


4nT Ju Zz 
. =. w(u, z/B) (32) 
where 
™* a3 (33) 
and 
we (2) (34) 


Cooper (1963) expressed the solution as: 
- - e”-(v*/ 
ad rcs l. y x dy 





= (35) 
ra L (u, v) 
with 
r x! 
ie (5) (36) 


The notations of Hantush and Jacob (1955) and Cooper (1963) are in- 
cluded here because tvpe curves, tables, and data analyses using both 
are found in the literature. Hantush and Jacob (1955) point out that as 
B approaches infinity, that is as leakage decreases, equation 32 ap- 
proaches the Theis equation (Equation 6). The L(u,v) of Cooper (1963) is 
called the leakance function of u and v. 


Hantush and Jacob (1954) noted that flow in a leaky confined aquifer is 
three dimensional, but if the hydraulic conductivity of the aquifer is 
sufficiently greater than that of the leaky confining bed, the flow may be 
assumed to be vertical in the confining bed and radial in the aquifer. 
This relation has been quantified by Hantush (1967a) for the condition 
b/B < 0.1. Assumption 5, that there is no change in storage of water in 
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Figure 13. Type curve of W(u, r / B) - L(u,v) as a function of 1 /u (From Loman, 1972). 
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the confining bed, was investigated by Neuman and Witherspoon 
(1969a). They concluded that this assumption would not affect the 
solution if 8 < 0.01, where 


6 = z [se 
4b NS, (37) 


Assumption 4, that there is no drawdown in water level in the source 
bed, was also examined by Neuman and Witherspoon (1969a). They 
indicated that drawdown in the source bed is justified when T, > 100T, 
where T, represents the transmissivity of the source bed and would have 
negligible effect on the drawdown in the pumped aquifer for short times; 
that is, when Tt/r’S < 1.6 B*(r/B)’. 


Figure 13 shows plots of dimensionless drawdown compared to dimen- 
sionless time from Reed (1980) using t/se notations of Hantush, Jacob, 
and Cooper (1963). 


ere 


Aquifer test data may be plotted in two ways. For the first method, 
measured drawdown in any one well is plotted versus t/r’; the data are 
matched to the solid line type curves in Figure 13. The data points are 
aligned with the solid-line type curves either on one of them or between 
them. Using the notation of Hantush and Jacob (1955), the parameters 
are then computed from the coordinates of the match points (t/r’,s) and 
[1/u, W(u,r/B)], and an interpolated value of r/B from the equations 


g=4T (242) (39) 


x’. T 
>» B (40) 


Using the notation of Cooper (1963), the parameters are computed from 
the coordinates of the match points (t/r’,s) and [1/u,L(u,v)], and an inter- 


polated value of v from equations 
re (2) (2iuv)) (41) 
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Figure 14. Type curve of the Bessel function K) x (x) as a function of x (After Reed, 1980). 
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2 
g=4T (He) (42) 


x’ v? 
> =» 4T (>) (43) 


This method was used by Cooper (1963) and the data and analysis of 
Cooper is cited by Lohman (1972). 


Cooper (1963) devised a second method as discussed by Bedinger and 
others (1988) by which drawdown measured at the same time but in 
different wells at different distances can be plotted versus t/r’ and 
matched to the dashed curves of Figure 13. The data are matched so as 
to align with the dashed line curves, either on one or between two of 
them. From the match point coordinates (s,t/r’) and [W(u,r/B), 1/u] and 
an interpolated value of v’/u, T and S are computed from equations 41 


and 42 and the remaining parameter from 


E es (ele) (44) 


Equilibrium Method relates to the fact that the zone v’/u 2 8 and 
W(u,r/B) 2 0.02 in the method of Hantush (1956) corresponds to steady 
state conditions. The drawdown in the steady state zone is given by 
Jacob (1946): 


ss (3%) K, (x) (45) 


where Ko(x) is the zero order modified Bessel function of the second kind 


x’ 
aoe (5) (46) 


Data for steady state conditions can be analyzed using the type curve in 
Figure 14. The drawdowns are plotted versus r and matched to the type 
curve. After choosing a convenient match point with coordinates (s,r) 
and [Ko(x),x], the parameters are computed from the equations: 


tT. (sh K, (x) (47) 


41 
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Figure 15. Section through a pumping well in a leaky confined aquifer with stor- 
age of water in confining beds (After Dawson and Istok, 1991). 
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£2 (48) 


Leaky Confining Bed With Storage 
Assumptions: 


1. Pumping well discharge is at a constant rate, Q. 


2. Pumping well is of infinitesimal diameter and fully penetrates 
the confined aquifer. 


3. Confined aquifer is overlain and underlain everywhere by 
leaky confining beds having uniform values of hydraulic 
conductivity, K’ and K”", thickness, b' and b", and storage 
coefficient, S' and S”. 


4. Leaky confining beds are overlain and underlain everywhere 
by infinite beds with constant hydraulic heads. 


5. Flow in the confined aquifer is two dimensional and radial in 
the horizontal plane; flow in the leaky confining beds is 
vertical. 


Hantush (1960) presented solutions for determining head in response to 
discharge from leaky confined aquifers where release of water from stor- 
age in the confining beds is taken into account. Release of water from 
storage in confining beds may be substantial in a number or geohydro- 
logic situations, such as where the confining beds are thick or where the 
upper confining bed contains a water table. Also, release of water from 
storage may be substantial for short durations (t less than both b' S’/10K’ 
and b" S"/10K") in many geohydrologic situations. Release of water from 
storage in confining beds commonly becomes less significant with time as 
steady state flow conditions are approached. A complete discussion of the 
Hantush (1960) methods for the geometry in Figure 15 and for other 
types of geometry is presented by Reed (1980). The solution of Hantush 
(1960) for short durations (t less than both b’ S'/10K' and b" S"/10K’) is: 


s- (22,) A (u, 6) (49) 
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Figure 16. Dimensionless drawdown, H a function of dimensionless time, 1/u, for a well fully penetrating 
a leaky confined aquifer wi SB danas deuntaen 1972). 
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where 

us rs (50) 
and 

p= 2 (Ee) « «(zeeth (51) 
and 

H(u,) =f” 2. erfc FG dy (52) 
and 


2 >. 
erfc(x) = Va f- e dy (53) 


Lohman (1972) points out that the versatility of equations 49 through 53 
is because they are the general solution for determining the drawdown 
distribution in all confined aquifers as discussed by Bedinger and others 
(1988), whether they are leaky or nonleaky. That is, 8 approaches zero 
as K’ and K" approach zero, and equation 48 becomes the Theis equation 9. 


beatinttion: 


The method can be applied by plotting drawdown versus t/r” and super- 
posing the data plot on the type curve plot of H(u,8) versus 1/u as shown 
in Figure 16. An example of the application of this method using data 
from an aquifer test is presented by Lohman (1972). 


Unconfined Conditions 


Conditions governing drawdown due to discharge from an unconfined 
aquifer differ markedly from those due to pumping from a nonleaky 
confined aquifer. Difficulties in deriving analytical solutions to the 
hydraulic head distribution in an unconfined aquifer result from the 
following characteristics: 
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Figure 17. Diagrammatic section through a pumping well in an unconfined 
aquifer (After Fetter, 1988). 
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1. Transmissivity varies in space and time as the water table is 
drawn down and the aquifer is dewatered. 


2. Water is derived from storage in an unconfined aquifer mainly 
at the free water surface and, to a smaller degree, from each 
discrete point within the aquifer. 


3. Vertical components of flow exist in the aquifer in response to 
withdrawal of water from a well. These components may be 
large and are greater near the pumping well and at early 
times. The diagrammatic section in Figure 17 is through a 
pumping well in an unconfined aquifer and shows conditions 
when the pumping well is pumped. If the water level in the 
pumping well is below the top of the screened interval, a 
seepage face will be present. 


The drawdown curve in an unconfined aquifer in response to an active 
pumping well follows a typical S-shaped curve. During early times of 
pumping activity, water level decline is rapid, and water is derived 
internally from the aquifer by expansion of the water and compaction of 
the aquifer; head response is similar to that of a confined aquifer. As 
pumping continues, head response lags that of a confined aquifer. This 
lag was attributed to slow drainage from the unsaturated zone by many 
early investigators. However, Cooley and Case (1973) concluded that the 
unsaturated zone has little effect on flow in the aquifer. Neuman (1972) 
attributed the lag to delayed response related to vertical components of 
flow in the aquifer as a function of the radial distance from the pumping 
well and of time. At later times, the drawdown once again appears to 
follow the Theis curve. 


Solution of Boulton and Neuman 


Boulton (1954b and 1963) introduced a mathematical solution to the 
head distribution in response to pumping an unconfined aquifer. 
Boulton’s solution derives the typical S-shaped curves of unconfined 
aquifers, but invokes the use of a semiempirical delay index that was not 
defined on a physical basis as discussed by Bedinger and others (1988). 
Neuman (1972 and 1975) presented a solution for unconfined aquifers 
based on well-defined physical properties of the aquifer. Neuman (1975) 
examined the physical basis for Boulton’s delay index (1/a) and deter- 
mined that, for fully penetrating pumping wells, Boulton’s solution 
yielded values of transmissivity, specific yield, and storage coefficient 
identical to those determined by Neuman. Neuman’s method for uncon- 
fined aquifers is discussed here. For further information on Boulton’s 
method, the reader is referred to Boulton (1954a and 1963). Application 
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of Boulton’s method to aquifer test data from unconfined aquifers is 
presented by Prickett (1965) and Lohman (1972). 


Assumptions: 


1. Pumping well of infinitesimal diameter discharges at a con- 
stant rate, Q. 


2. Pumping well and observation well are open throughout the 
thickness of the unconfined aquifer. 


3. Unconfined aquifer is areally extensive, homogeneous, and 
isotropic with vertical hydraulic conductivity, K,, and horizon- 
tal hydraulic conductivity, K,. 
Solution: 


The solution of Neuman (1973) for the condition in which the pumping 
well and the observation well are perforated throughout the saturated 
section of the aquifer is given by: 


s(r,t) = 22 J, 47% (ve™) 





[uy (y) p> u,(y)] dy (54) 
where 
u(y) = Soe ee ean. its 
and 
u(y) = £i2Bxp [ot,B (y?-y9) 1) tan (y_) “ 





{y?- (1+) y5-[(y?-75)?/0)} y, 
and the terms Y and yy, are the roots of the equations 


oy, sinh (y,) - (y?-y%) cosh (y,) = 0 


where y, < y? (57) 
shows relationship of quantities 
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oy, sin (y,) + (y?+y3) cos (y,) = 0 
@ (2n-) (#/2) < y, < ae, nz1 (58) 
shows relationship of quantities 


Equations 54 through 56 are expressed in terms of three independent 
dimensionless parameters, 8, tz, and s. Neuman (1975) decreased the 
number of independent dimensionless parameters by considering the 
case in which s=S/Sy approaches zero, that is, in which S is much less 
than Sy. The results are two asymptotic families of type curves referred 
to as type A and type B curves (Figure 18). Neuman (1975) listed nu- 
merical values for the curves. 


The curves lying to the left of the values of 8 in Figure 18 are called type 
A curves and correspond to the top scale expressed in terms of ts. The 
curves lying to the right of the values of 8 in Figure 18 are called type B 
curves and correspond to the bottom scale expressed in terms of ty. The 
two sets of curves are asymptotic to Theis curves. Type A curves are 
intended for use with early drawdown data and type B curves with late 
drawdown data. 


sities 


Neuman (1975) described application of his solution for aquifer charac- 
teristics by two methods: Using logarithmic plots of aquifer test data 
and type curves, and using semilogarithmic plots of aquifer test data. 


The logarithmic method as described by Neuman (1975) follows. Late- 
time drawdown, s, is plotted for the observation well on logarithmic 
tracing paper against values of time, t. This data plot is overlain on the 
type B curves; while keeping the vertical and horizontal axes of both 
graphs parallel, as much of the late time drawdown data is matched to a 
particular curve as possible and a match point is selected. The value of 
B of the type curve matched is noted and the coordinates of s, sp, and t, 
ty of the match point are recorded. The transmissivity is calculated from 





and the specific yield from 
Tt 
sS.«s 
7 (60) 


Next, the process is repeated by overlaying the early time drawdown 
data on the type Acurves. The value of 8 corresponding to the type 
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Figure 18. Type curves for fully penetrating wells in an unconfined aquifer (After Neuman, 1975). 
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curve must be the same as that obtained earlier from the B curves. The 
coordinates of the new match point s, sp, and t, t, are recorded. The 
transmissivity is again calculated from equation 59. Its value should be 
approximately equal to the previously calculated value from the late-time 
drawdown data as discussed by Bedinger and others (1988). The storage 
coefficient is calculated from: 


Tt 


o° ie. (61) 


The horizontal hydraulic conductivity, K,, is calculated from the value of 
B according to: 


Kz (62) 
The degree of anisotropy, Kp is calculated from: 
2 
K, = 22 (63) 
The vertical hydraulic conductivity, K,, is calculated from: 
K, = K,K, (64) 


Utilization of Confined Aquifer Methods to Unconfined Aquifers 


The methods of Theis (1935) ar.d Theim (1906), and other methods, 
though applicable to confined aquifers, may also be applied to unconfined 
aquifers where the drawdown is small in relation to the thickness of the 
aquifer (Jacob, 1950). Corrections in drawdown need to be made when 
the drawdown is a significant fraction of the aquifer thickness. Such 
corrections are usually called thin-aquifer corrections. These methods 
rely on the Dupuit-Forcheimer assumptions and are not valid for early 
time when vertical flow components are substantial. The Dupuit- 
Forcheimer assumptions state that, within the cone of depression of a 
pumping well, the head is constant throughout any vertical line through 
the aquifer and is, therefore, represented by the water table as discussed 
by Kruseman and Ridder (1991). Actually, this is true only in a confined 
aquifer having uniform hydraulic conductivity and a fully penetrating 
pumping well. Jacob (1963a) stated that where the drawdown needs to 
be replaced by s', the drawdown that would occur in an equivalent con- 


fined aquifer would be represented by 
s? 
a'=s4- (>5) (65) 


Jacob (1963a) presented a correction for the coefficient of storage where 
the drawdown is a substantial fraction of the original saturated 
thickness, 
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g's (2-8) 3 (66) 


Neuman (1975) recommended that the use of confined aquifer methods 
for unconfined aquifers be restricted to late time data after the effect of 


delayed gravity response. 


Where vertical components of head may be substantial, paired observa- 
tion wells that partially penetrate the aquifer may be used in lieu of 
fully penetrating observation wells. One well of the pair is screened at 
the bottom of the aquifer and the other is screened just below the water 
table. The water levels in the paired wells are averaged and used in the 
confined aquifer method along with the thin aquifer correction, as neces- 
sary (Lohman, 1972). 


Estimating Stream Depletion by Pumping Wells 


The correlation between stresses imposed by pumping wells and the 
resultant depletion of stream flows has been identified by numerous 
investigators (Glover and Balmer, 1960; Theis and Conover, 1963; 
Hantush, 1964). This correlation is usually shown by charts and equa- 
tions as discussed by Jenkins (1970). The techniques shown in this 
section are mainly derived from the work of Jenkins (1970) who provided 
easy to follow tools such as curves, tables, and sample computations. 
The symbols that are employed in this section are defined below: 


transmissivity [L’*/T} 

specific storage of the aquifer, dimensionless 

time, during pumping period, since pumping began [T] 
total time of pumping 

time after pumping stops [T] 

net steady pumping rate [L’/T] 

rate of depletion of the stream [L’/T] 

net volume pumped during time t [L’] 

net volume pumped [L’] 

volume of stream depletion during time t, tp + ti [L’] 
perpendicular distance from pumped well to stream [L] 
stream depletion factor [T]. 


Jenkins (1970) defines the stream depletion factor to be the time coordi- 
nate where the volume of stream depletion is equal to 28 percent of Qt 
on a curve relating v to t, and if sdf = a’S/T. In a complex system, it can 
be considered to be an effective value of a®°S/T. Jenkins (1970) further 
states that the value of the sdf at any location in the system depends 


upon the integrated effects of the following: irregular impermeable 


2” ‘or 
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tion, distance from the stream, and specific hydraulic connection be- 
tween the stream and the aquifer. It should be noted that the curves 
and tables used for calculating depletions in this section are dimension- 
less and can be used with any units as long as they are consistent. 
Jenkins (1970) states that the assumptions used in the analysis of calcu- 
lating stream depletion from pumping wells are as follows: 


1. T does not change with time. 

2. The temperature of the stream is assumed to be constant and 
the same as the temperature of the aquifer. 

3. The aquifer is isotropic, homogenous, and semi-infinite in 
areal extent. 

4. The stream that forms the boundary is straight and fully pen- 
etrates the aquifer. 

5. Water is released instantaneously from storage. 

6. The well is open to the full saturated thickness of the aquifer. 

7. The pumping rate is steady during any period of pumping. 


Curves A and B in Figure 19 apply during the period of steady state 
pumping as discussed by Jenkins (1970). Curve A defines the correla- 
tion between the dimensionless term, t/sdf, and the rate of stream deple- 
tion, q, at time t, and is shown as a ratio to the pumping rate Q. Curve 
B defines the correlation between t/sdf and the volume of the stream 
depletion, v, during time t, and is defined as a ratio to the volume 
pumped, Qt. The curves 1-q/Q and 1-v/Qt are defined to better interpret 
values of q/Q and v/Qt when the ratios surpass 0.5. The coordinates of 
curves A and B are tabulated in Table 2. The curves A and B that are 
tabulated in this section are after Jenkins (1970). It should be noted 
that the precision is only to two significant places, which is considered to 
be appropriate for this type of analysis. 


Sample Problem 


To explain the application of the curves and table, a sample problem is 
defined and solved using the method outlined in this section. The prob- 
lem is typical of what may be encountered in the field or as a proposed 
activity. It is assumed that the data used in the examples is usually 

available during the field study phases of most water resource projects. 


The problem is a pumping well that is pumped at 2.0 acre feet per day 
and is located 1.58 miles from a stream. The question is: How long can 
the pumping continue before the stream depletion reaches 0.14 acre feet 
per day, and what is the total stream depletion for the period of pumping? 







































































Ppeeecres 


























i 




















Con 


i. 



































10 























t/sdf 


1.0 

















(ma eeer 
[eens 






































> FP 
Curve A (@/Q) 


[= aes 






























































1.0 





0.1 
0.01 


1/A GNV Of) 





























0.1 








(After Jenkins, 1970). 


























0.001 VL 
0.001 





Figure 19. Curves to interpret rate and volume of stream depletion 





Table 2. Values of q/Q, Q, v/Qt, and v/Qsdf relating selected values of t/sdf 


Values of g/Q, v/Qt, and v/Qsdf corresponding -o selected values of t/sdf 
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0.14 ac-ft/d 

2.0 ac-fi/d 

1.58 miles 

106 gal/d/ft 

30 days 

a’S/T = a’(T/S) = 
(1.58mi)’ (5,280 ft/mi)’ / 10° gal/d/ft) (1 ft"/7.48 gal) = 
520 days 





From the data given, the ratio of the rate of stream depletion to the rate 
of pumping is q/Q = (0.14ac-ft/d(2.0 ac-ft/d) = 0.07 


From curve A (Figure 19) t/sdf = 0.15 


Substitute the value under Data for sdf, and t = (0.15) (520 days) = 78 
days 


The total time the well can be pumped is 78 days 
When t/sdf = 0.15 then from curve B (Figure 19) v/Qt = 0.02 


Substitute the values for Q and t, and the volume of the stream depletion 
during this time period is v = (0.02) (2ac-ft/d) (78 days) = 3.1 ac-ft 


This shows that during the 78-day purnping period, 3.1 ac-ft of water can 
be attributed as stream depletion. 


It should be noted that variation from idealized stream conditions may 
cause actual stream depletions to be either more or less than the values 
interpreted from the method discussed in this section. Fluctuations in 
water temperature will cause variations in stream particularly 
by large-capacity wells near stream lengths. As discussed by Moore and 
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Jenkins (1966), if large-capacity wells are located close to a stream 
length and streambed permeability is low compared to aquifer perme- 
ability, the water table may be drawn down below the bottom of the 
streambed. The methods discussed in this section are not appropriate 
for streambed permeability, area of the streambed, temperature of the 
water, and stage of the stream. 


The mathemw‘ic#) basis for the curve development and table presented 
in this section is beyond the scope of this guide. If the reader is inter- 
ested in a more detailed discussion of the mathematical curve and table 
development, they are referred to the work of Glover (1954), Jenkins 
(1968a), Theis (1941), and Theis and Conover (1963). 
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Pumping Well and Flow System 
Characteristics 


Pumping Well Conditions 


The aquifer test methods discussed in the previous sections are, for the 
most part, based on simple geometric conditions and constant discharge 
or head change in the pumping well. For example, in the methods for 
confined aquifers, the pumping well was assumed to fully penetrate the 
aquifer and result in radial flow and vertically uniform heads in the 
aquifer. In the methods for unconfined aquifers, the pumping well and 
the observation well were assumed to be fully penetrating, which simpli- 
fies the analytical solution and application of the method to the problem. 
Storage effects in the pumping well were assumed to be negligible, ex- 
cept in the slug test methods. Variations from simple geometric condi- 
tions of the pumping well and variable discharge may cause anomalous 
hydraulic conditions in the flow field as discussed by Bedinger and 
others (1988). Disregard of pumping well geometry, that is, such factors 
as length of screened interval and depth of screen, or in some caves, uf 
observation well geometry, may invalidate an aquifer test method. In 
this section, some principal methods are discussed for accounting for 
partially penetrating wells, variations in the discharge rate of the pump- 
ing well, constant drawdown, storage in the pumping well, inertial 
effects of water in the pumping well and in the aquifer, multiple aqui- 
fers, fractured media, anisotropic media, and image wells. 


Partially Penetrating Wells 


Partially penetrating pumping wells cause vertical components of flow 
that greatly complicate the analytical solution to the hydraulic head 
distribution in the aquifer and the application of the solution to aquifer 
tests. For these reasons, methods of confined aquifer test analyses 
treated previously assume fully penetrating pumping wells and radial 
flow. The analytical method of Neuman (1972 and 1975) for unconfined 
aquifers presented in the section on Hydraulic Test Methods for Aquifers 
is based on the assumptions of completely penetrating pumping and 
observation wells. These assumptions made it possible to simplify the 
application of the method using a single family of type curves. 


The vertical components of head caused by partial penetration of the 

pumping well need to be considered for (K,/K,)” r/b < 1.5 (Reed, 1980). 
Thus, in a homogeneous, isotropic confined aquifer, where K, = Kr, the 
effects of a partially penetrating pumping well are negligible beyond a 
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distance of about 1.5 times the aquifer thickness as discussed by 
Bedinger and others (1988). In an aquifer having radial-vertical anisot- 
ropy, where K, > K,, vertical flow components are of concern for a 
greater distance from the pumping well. Analytical solutions have been 
presented for anisotropic confined and unconfined aquifers and methods 
have been developed for their application to aquifer tests. These meth- 
ods are discussed in the section on Anisotropic Aquifer Materials. 


Variably Discharging Wells 


Stallman (1962) and Moench (1971) presented methods of analysis of 
drawdown in response to an arbitrary discharge function. These meth- 
ods simulate pumpage as a sequence of constant rate step changes in 
discharge. The methods utilize the principle of superposition in con- 
struction type curves by summing the effects of successive changes in 
discharge. The type curves may be derived for pumping wells discharg- 
ing from extensive, leaky, and nonleaky confined aquifers, or any situa- 
tion where the response to a unit stress is known. 


Recognizing that the uncontrolled discharge from a pumping well com- 
monly decreases with time during the early period of pumping, type 
curves have been described for drawdown in response to decreasing 
discharge functions that can be expressed mathematically. Hantush 
(1964b) developed drawdown formulas for three types of decrease in 
pumping well discharge including an exponentially decreasing discharge 
and a hyperbolically decreasing discharge for extensive, uniform con- 
fined aquifers. Methods for leaky, sloping leaky, and nonleaky confined 
aquifers also were presented by Hantush (1964b). 


Abu-Zied and Scott (1963) presented a general solution for drawdown in 
an extensive confined aquifer in which the discharge of the pumping 
well decreases at an exponential rate. Aror 9nd Scott (1965) proposed 
an approximate method of determining tri 1:1) 3sivity ard storage from 
an aquifer test in which discharge decreas. » th time during the early 
part of the test. 


Lai and Su (1974) presented methods for determining the drawdown in a 
homogeneous, isotropic, nonleaky, confined aquifer, taking into account 
storage in the pumping well in response to exponentially and linearly 
decreasing discharge. Lai and Su (1974) also presented a method for 
determining drawdown in a homogeneous, isotropic, leaky confined 
aquifer, taking into account storage in the pumping well and various 
discharge rates. 
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Constant Drawdown Conditions 


Methods to determine the hydraulic head distribution around a pumping 
well in a confined aquifer with near constant drawdown are presented 
by Jacob and Lohman (1952), Hantush (1964b), and Rushton and 
Rathod (1980). Such conditions are most commonly achieved by shut- 
ting in a flowing well long enough for the head to fully recover, then 
opening the well. The solutions of Jacob and Lohman (1952), and 
Hantush (1964a) apply to areally extensive, nonieaky confined aquifers. 
Rushton and Rathod (1980) used a numerical model to analyze aquifer 
test data. When using the method of Jacob and Lohman (1952), mea- 
surements are made of the decreasing rate of flow after the pumping 
well is opened. Application of the method by type curve and straight- 
line techniques is described by Lohman (1972). Hantush (1959b) pre- 
sented two methods for determining constant drawdown in a leaky 
confined aquifer without storage in the confining beds. One method is for 
the discharge of the pumping well; the other is for the drawdown in the 
aquifer. Reed (1980) presented a computer program for calculating 
function values for Hantush’s (1960) methods. The method of Hantush 
(1964a) uses measurements of head in the flowing pumping well and in 
an observation well to determine diffusivity (T/S). 


Storage and Inertial Influence 


The effect of storage and inertial effects of head in the pumping well and 
the aquifer were examined by Bredehoeft and others (1966). For con- 
tinuous pumping under ordinary conditions of pumping from production 
wells in transmissive aquifers, the effects of storage in a production well 
become negligible in a short time. However, the effects of storage could 
be significant in pumping wells of large diameter drawing water from 
aquifers having minimal transmissivity as discussed by Bedinger and 
others (1988). The effect of slug well storage is commonly significant in 
slug tests and the slug test methods presented in an earlier section of 
this guide account for storage in the slug well. Most aquifer test meth- 
ods do not consider the effect of storage within the pumping well; hence, 
the stated assumption that the pumping well is of infinitesimal diam- 
eter. According to Papadopulos and Cooper (1967), the pumping well 
storage may be neglected if t > 2.5 x 10°r,’/T, where r; is the radius of 
the well casing in the interval in which the water level declines. 


Papadopulos (1967b) presented a solution for determining the drawdown 
in and around a pumping well of finite diameter taking into consider- 
ation the effect of water stored in the wellbore. Papadopulos (1967b) 
presented tables and type curves and discussed application of type curve 
techniques for solution(s) to the problem. Tables and type curves for 
application of the method are presented in Reed (1980). 
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Inertial effects in a well are a function of the well and the aquifer. 
Force-free oscillation occurs in underdamped wells following events such 
as earthquakes or sudden imposition of a head change. Bredehoeft and 
others (1966) presented examples in which the column of water in 
underdamped wells oscillates for a few seconds after a sudden com- 
mencement of continuous pumping. Inertial effects to continuous pump- 
ing are probably not significant in most aquifer tests. 


Van der Kamp (1976) and Kipp (1985) presented methods for determin- 
ing the transmissivity of an aquifer from inertially induced oscillation in 
a pumping well, a response that may occur in conjunction with ex- 
tremely transmissive aquifers. Van der Kamp (1976) suggested a tech- 
nique for inducing oscillations in a well by a procedure used in some slug 
tests; that is, by sudden removal of a closed cylinder of known volume for 
the well. Kipp (1985) presented the complete method from the 
noninertially induced slug well response to the freely oscillating slug 
well. The method of Kipp (1985) is a useful extension to conventional 
slug test methods. Slug test methods are suitable for damped slug wells, 
those in which force-free oscillations are negligible as is common in 
aquifers with minimal to everage transmissivity. 


Flow System Characteristics 


Flow system characteristics for the aquifer test methods discussed in the 
previous sections were based on simple geohydrologic characteristics. 
Aquifers were assumed to be homogeneous and isotropic and of infinite 
areal extent. In this section, methods are introduced which deal with 
multiple, fractured, and anisotropic aquifers, and with aquifers of finite 
areal extent bounded by impermeable and constant head boundaries as 
discussed by Bedinger and others (1988). 


Multiple Aquifers 


Tests of multiple aquifers, that is, two or more aquifers separated by a 
leaky confining bed or penetrated by a pumping well, require special 
methods for analysis. Bennett and Patten (1962) devised a method for 
testing a multiaquifer system by a procedure using downhole metering 
and constant drawdown. Extending his work with leaky aquifers, 
Hantush (1967b) presented a solution for determining drawdown distri- 
bution in two aquifers separated by a leaky confining bed, in which 
storage is neglected, in response to discharge from one or both of the 
aquifers. Neuman (1972) provides a solution for drawdown in leaky 


confining beds above and below an aquifer being pumped. Neuman and 
Witherspoon (1969a) developed an analytical solution for the flow in a 
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leaky confined system of two aquifers separated by a leaky confining bed 
with storage. One of the aquifers is discharged through a fully penetrating 
well. Javandel and Witherspoon (1969) presented a finite element method 
of analyzing anisotropic multiaquifer systems. 


Fracture Flow 


Models that have been developed for flow in fractured rock include those 
based on the assumptions that flow is in a single fracture composed of 
parallel plates, flow is in a network of intersecting frectures, and flow is in a 
double porosity medium consisting of blocks containing intergranular poros- 
ity and permeability, the blocks being separated by « network of intersect- 
ing fractures sufficiently extensive to be considered a continuum. A review 
of methods of treating fractured media is presented in Gringarten (1982). 


Solution for flow in single finite fractures in a porous medium is presented 
by Gringarten and Ramey (1974). Barenblatt and others (1960) presented a 
method for solving the double porosity model. This model is based on the 
assumptions that storage of water in the fractures is negligible compared to 
storage in the pores of the blocks, and flow of water is primarily in the 
fractures. Boulton and Streltsova (1977) presented a solution for a system 
composed of porous layers separated by fractured layers that are horizontal. 
Moench (1984) developed type curves for a double porosity model with a 
fracture skin that may be present at the fracture block interfaces as a result 
of mineral deposition or alteraticn. 


Anisotropic Aquifer Materials 


Most of t'ne aquifer test methods discussed in the section on H)draulic Test 
Methods for Aquifers are based on the assumption that the aquifer is homo- 
geneous and isotropic. Natural materials are neither homogeneous nor 
isotrozic, but aquifer test methods based on this assumption aie widely 
applied and useful. Aquifer test methods and associated procedures that 
have been devised to evaluate anisotropy of natural media will be intro- 
duced in this section. Vertical anisotropy is common in stratified sediments. 
Anisotropy often is characteristic of natural formations. Hantush (1966a 
and 1966b) presented methods for determining flow in homogeneous, aniso- 
tropic media, but did not provide procedures for applying the methods. 
Methods described in the literature for treating anisotropy are limited to 
the situations of either horizontal anisotropy or horizontal and vertical 
anisotropy. 

Solutions to the head distribution in a homogeneous confined aquifer with 


radial vertical anisotropy in response to constant discharge of a partially 
penetrating well are presented by Hantush (1961a). 
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Figure 20. Section showing drawdown and flow paths near a pumping well in an 
ideal nonleaky confined aquifer. Solid lines are drawdown and flow 
lines for a pumping well screen in bottom of aquifer; dashed lines are 
drawdown lines for a pumping well screened the full aquifer thickness 
(After Weeks, 1969). 
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The solutions of Hantush (1961a) were applied by Weeks (1964 and 
1969), who presented methods to determine the ratio of horizontal to 
vertical hydraulic conductivity. The analyses are made by comparing 
measure drawdowns in the piezometers to those predicted if the pump- 
ing well fully penetrates the aquifer. The differences in the measured 
partially penetrating pumping well at which these differences would 
occur in an isotropic aquifer are determined from an equation. The 
permeability ratio is computed as the square of the ratio of the actual 
distances to the computed distances. Weeks (1969) applied graphical 
methods to the solution of vertical and horizontal hydraulic conductivity 
and prese ited tables of values of the dimensionless drawdown correction 
factor (Figure 20). Weeks (1969) also discussed conditions for which his 
method is applicable to unconfined aquifers. 


Papadopulos (1965) presented a method for determination of horizontal 
plane anisotropy in an areally extensive, homogeneous, confined aquifer. 
Papadopulos (1965) introduced a graphical method for solution of the 
components of the transmissivity tensor from aquifer test data using a 
minimum of three observation wells. Hantush and Thomas (1966) pre- 
sented a graphical method of determining horizontal anisotropy in con- 
fined aquifers from the elliptical shape of the cone of drawdown. 


Neuman (1975) presented a solution for the drawdown in piezometers in 
response to discharge from a partially penetrating pumping well in an 
unconfined aquifer having radial and vertical anisotropy. Because of the 
large number of variables involved, Neuman (1975) offered to provide a 
computer program from which the user could prepare type curves for 


specific cases. 
Image Well Method 


Each of the aquifer test methods of aquifer tests discussed previously in 
this guide is based on the assumption that the aquifer is of infinite areal 
extent. It is recognized that such conditions do not exist. Effects of 
limitations in areal extent of aquifers by impermeable boundaries or by 
source boundaries, such as hydraulically connected streams, may pre- 
clude the direct application of an aquifer test method. The method of 
images provides a tool by which a solution to the problem of exterior 
boundaries can be devised as discussed by Bedinger and others (1988). 
This method uses the substitution of a hydraulic boundary for the physi- 
cal feature. 


Consider first an aquifer bounded by a perennial stream in which the 
head is independent of the pumping well; that is, there is no drawdown 
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in the stream and the stream functions as a fully penetrating, constant 
head boundary to the aquifer (Figure 21A). An image system that satis- 
fies the foregoing boundary condition is shown in Figure 21B; that is, an 
imaginary recharging well located on the opposite side of and the same 
distance from the stream as the real pumping well. Both wells are on a 
line perpendicular to the stream. The imaginary recharge well operates 
simultaneously with the real pumping well and recharges water to the 
system at the same rate the real well discharges. The resultant draw- 
down at any point in the system is the algebraic sum of the drawdown 
ES eS COS SES GS a Oa SS Se se ey 
nary wells. 


Next, consider an aquifer bounded by confining material (Figure 22A). 
The hydraulic boundary condition imposed by the confining material is 
that there is no flow across the material. The image well condition that 
duplicates this physical condition by hydraulic analogy is shown in 
Figure 22B. An imaginary pumping well has been placed at the same 
distance from the line of zero flow as the real well. The wells are on the 
opposite sides of and on a line perpendicular to the line of zero flow as 
discussed by Bedinger and others (1988). As in the case of the flow 
system with the recharging image well, the resultant drawdown at any 
point in the system is the algebraic sum of the changes in head caused 
by the real and imaginary well. 


The theory of images may be applied to any combination of straight-line 
constant head and impermeable boundaries. A number of combinations 
are discussed by Ferris and others (1962). Because the drawdown in an 
observation well, 8s, in a system bounded by a line source or imperme- 
able boundary, is the algebraic sum of the components of drawdown by 
the pumping well, 8p, and by the image well, sj, the hydraulic head 
distribution in the aquifer can be analyzed by superposing the solutions, 
by an appropriate aquifer test method for the flow system, for the real 
and image wells. For example, if the flow system is a confined nonleaky 
aquifer, it would be appropriate to apply the method of Theis in 
equation 4. The Theis equation for an aquifer Sounded by line source or 
impermeable boundary where 

5° 8,2 8; (67) 


becomes (Stallman, 1963) 


s, = =f. (W(u) p + W(u) i} 


= (68) 
(ks > W(u) 
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Figure 22. Idealized sections of a pumping well in a semi-infinite aquifer 
bounded by an imperineable formation and of the equive’ ~nt nvdralic 
system in an infinite aquifer (After Ferris and others, 1962). 
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and 

u,* eS (69) 
and 

u,* zi (70) 


Type curves can bi constructed for a specific observation well or a family 
of type curves can be drawn for different ratios of rj = rp. Such a family 
of type curves is presented by Stallman (1963) and Lohman (1972), who 
discuss application of the method. The type curves can be used to ana- 
lyze the drawdown data for hydraulic properties of the aquifer in a 
system where a boundary is known to occur or to locate the position of a 
hidden boundary that is indicated by the drawdown data from an aqui- 
fer test (Morris and others, 1959; Moulder, 1963). Boundaries in nature 
may be neither absolutely impermeable nor constant head. For ex- 
ample, streams generally do not fully penetrate the aquifer and 
streambed materials :~ay limit the rate of water movement from the 
stream to the aquifer. 


Methods to determine an effective distance from a pumped well to a 
stream boundary include a type curve method suggested by Kazman 
(1946) that is implicit in the type curves of Stallman (1963) and Lohman 
(1972), and in the graphical extrapolation of the drawdown to zero from 
the weer level in a line of observation wells perpendicular to the river 
from the pumped well of Rorabaugh (1956). 
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Summary 


Geohydrology is an important part of all ground water resource projects 
with respect to analysis and design. Water resource projects include the 
modeling, planning, analysis, and interpretation of information on the 
subsurface environment of ground water. Critical elements of the data 
collection phases of a ground water resource project may be closely 
associated with geohydrologic testing. 


The expanding scientific literature on ground water hydraulics and 
hydrology should be read and evaluated on a continuing basis. The 
information contained in this guide will need to be supplemented as 
more updated methodologies and techniques become available. 
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Appendix A 
Summary of Hydraulic Test Methods 


This summary includes methods for hydraulic testing classified by aqui- 
fer condition, pumping well characteristics, recharge and discharge 
function, and boundary conditions. The summary is divided into three 
parts: Confined Aquifer, Unconfined Aquifer, and Other Conditions. 


I. Confined Aquifer 
A. Nonleaky confined aquifer 


Methods included here are for radial flow in a nonleaky, 
porous, homogeneous, and isotropic medium of infinite 
areal extent. Change in water stored is instantaneous and 
proportional to the change in head. The aquifer is confined 
above and below by impermeable beds. The water level is 
above the top of the aquifer. 


1. Constant flux. 


Theim (1906) - Asymptotic (pseudosteady) solution 
Theis (1935) - Negligible storage in pumped well 


Cooper and Jacob (1946) - Asymptotic (logarithmic) 
approximation to well function of Theis (1935) with 


increasing time and decreasing radial distance. 


Stallman (1963) - Aquifer bounded on one side by a 
straight boundary (either constant head or no flow.) 


2. Constant drawdown. 


Jacob and Lohman (1952) - Step change in water 
level at the pumping well. Discharge of the pumping 
well as a function of time. Commonly used for shut- 
in flowing pumping wells. 








; 
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Instantaneous head change (slug tests). 


A rapid change in water level is induced in the slug 
well by various methods, such as injection, bailing, or 
pressurization. Inertial effects are assumed to be 
negligible. Inertially induced oscillatory fluctuations 
and applications to slug tests are treated by Krauss 
(1974), Van der Kamp (1976), Shinohara and Ramey 
(1979), and Kipp (1985). 


Hvorslev (1951) - Applies differential equation for 
permeameters to head change in an aquifer. Cases 
involving both radial and vertical flow are treated by 
shape factors for different flow geometries. 


Skibitzke (1958) - A method for determining the 
water level in a well after it has been bailed. Bailed 
well is assumed to be a fully penetrating line source 
rather than a well of finite diameter. 


Ferris and Knowles (1963) - Change in water level is 
caused by a sudden injection of water. Injection well 
is assumed to be a fully penetrating line source 
rather than a well of finite diameter. 


Cooper and others (1963) - Derives equation, pre- 
sents curves and a table of functions, presents a 
method of determining transmissivity and storage 
coefficient taking into account well storage, and 
discusses the relation to the solution of Ferris and 
Knowles (1963). Can be applied to fractured rock 
(Wang and others, 1977) if fracture openings do not 
change with pressure and there is negligible drain- 
age form the matrix into the fractures. 


Papadopulos and others (1973) - Presents additional 
function values and curves for the method of Cooper 
and others (1967). 


Bredehoeft and Papadopulos (1980) - Discusses 
testing formations with minimal permeability by 
pressurizing a shut-in well. For a certain range of 
parameter values, the method of Cooper and others 
(1967) indicates only the product of transmissivity 
and storage. 
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Neuzil (1982) - Discusses changes in procedure and 
equipment for the method of Bredehoeft and 
Papadopuos (1980). 


Barker and Black (1983) - Considers an aquifer with 
uniform horizontal fissures, horizontal flow in the 

fissures, vertical flow in the matrix, and storage in 

both. Numerical inversions are used for analysis of 
errors resulting from matching slug test data form a 
fissured aquifer by the method of Cooper and others 
(1967), which was developed for a homogeneous aqui- 
fer. It was concluded that transmissivity will always 
be overestimated, although not likely by more than a 
factor of 3. Storage coefficient, however, could be 

several orders of magnitude either larger or smaller. 


Dougherty and Babu (1984) - Application of slug tests 
to fractured porous aquifers. 


Variable discharge. 


Flux of pumping well is not constant, but varies with 
time. 


Werner (1946) - Flux is a linear function of time. 


Stallman (1962) - Continuously varying discharge is 


approximated by step changes. Function curves for 
drawdown are sums of well function (Theis, 1935) 


weighted by the change in discharge. 

Abu-Zied and Scott (1963) - Flux exponentially 
changes with time. 

Abu-Zied and others (1964) - Treats special cases that 
simplify the method of Abu-Zied and Scott (1963). 


Aron and Scott (1965) - Superposes the log asymptote 
to the solution of Theis (1935) weighted by the change 


in discharge. 


Sternberg (1968) - Graphical summation based on the 
log approximation to the well function and a multiple- 
step approximation of the well discharge. 
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Moench (1971) - Convolution integral applied to a 
general discharge function and the method of Theis 
(1935). Representation of discharge as a step curve 
and evaluation of the integral by summation. 


Lai and others (1973) - Includes the effect of storage 
in the well. Presents a general solution as a convolu- 
tion integral. Presents solutions for exponential and 
linear discharge as unevaluated integrals of complex 
functions. No tables of values. Presents three type 

curves for linear decreasing flux of drawdown in the 


pumping well. 
Multiple Aquifers 


Papadopulos (1966) - Two nonleaky aquifers, with 
different hydraulic properties, separated by a confin- 
ing layer. Constant discharge from a well open to 


B. Nonleakv, fractured, confined aquifers. 


Fractures, rather than the medium, transmit most of the 
fluid, especially in the vicinity of the pumping well. 
Gringarten (1982) - Review articles discuss several aspects 
of flow to wells through fractured media. 


1. 


Extensive fractures. 


An extensive network of fractures, sufficiently dense 
and uniform as to be considered a continuum. 


Barenblatt and others (1960) - Double porosity model. 
Fractures in a porous medium. All storage in the 
pores. Flow from the medium into the fractures is 
proportional to the difference in head. Solution for 
head in the fractures. 

Warren and Root (1963) - Solutions to the conditions 


of Barenblatt and others (1960) that are applicable for 
long durations and for infinite and circular aquifers. 


Boulton and Streltsova (1977) - Radial flow in 
fractures that are separated by uniform layers of 
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porous medium. Only vertical flow within the layers. 
Storage in fractures and layers. Methods for frac- 
tures and layers. 


Dougherty and Babu (1984) - Analysis of slug tests in 
single and double-porosity aquifers by partially and 
fully penetrating wells with and without skin effect. 


Moench (1984) - Double-porosity model with fracture 
skin at fracture-block interfaces. 


Hsieh and others (1985) - Determination of three- 
dimensional, hydraulic-conductivity tensor in aniso- 
tropic fractured media. 


2. Single fracture. 
A single fracture centered about the pumping well. 
Gringarten and Ramey (1974) - Horizontal fracture. 
Gringarten and others (1974) - Vertical fracture 


Nonleaky confined aquifer with radial flow and horizontal 
anisotropy. 


Permeability and hydraulic conductivity are second order 
tensors in the horizontal plane. In two directions, the axes 
of the ellipse, flow, and hydraulic gradient are colinear. In 
other directions, flow and gradient are not parallel. 


Papadopulos (1965) - A minimum of ‘hree observation wells 
at different directions from the pur ~ »¢ well are needed to 
determine the principal comp...» and orientation of the 
transmissivity tensor. 


Hantush (1966a) - Drawdowns measured in three lines of 
observation wells are needed. A line is one or more obser- 
vation wells all in the same direction from the pumping 

well. If the principal directions are known, then two lines 


will permit analysis. 


Hantush (1966b) - By a simple transformation of coordi- 


nates, the boundary-value problem describing the flow in a 
homogeneous media may be transformed to an equivalent 
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homogeneous and isotropic aquifer. Methods are discussed 
for leaky confined aquifers, for complete and partial pen- 
etration, and for decreasing discharge. 


Hantush and Thomas (1966) - Distribution of observation 
wells is such that the elliptical shape of an equal drawdown 
(or residual drawdown for recovery) contour can be defined. 


Leaky and confined aquifer with radial flow and isotropic 
and homogeneous porous media. 


Vertical flow in uniform confining beds. Change in water 
stored is instantaneous with and proportional to change in 
head. Aquifer is confined above and below. Water level is 
above the top of the aquifer. Change in flow between aqui- 
fer and confining beds is proportional to drawdown. 


1. Constant flux 


Jacob (1946) - Solutions for steady flow in an exten- 
sive aquifer and nonsteady flow to a well at the 
center of a circular aquifer with no drawdown at the 
outer boundary. 


Hantush and Jacob (1955) - Solution for nonsteady 
flow in an extensive aquifer. 


Hantush (1956) - Graphical methods are applied to 
determine parameters. 


Hantush (1960) - Solutions that apply at either short 
or long durations and that include the effect of stor- 
age in confining beds. The three combinations of 
zero drawdown of zero-flow boundaries above and 
below the system are presented. 


Moench (1985) - Solution to transient flow to a well 
in an aquifer accounting for storage in the well and 
in the confining beds. 


Multiple aquifers. 

Hantush (1967b) - Two aquifers separated by a con- 
fining bed. No storage in the confining bed. Radial 
flow in the aquifers and vertical leakage through the 
confining bed is considered. Constant discharge from 
one aquifer. 


Neuman and Witherspoon (1969a) - Radial flow in 
two aquifers separated by a confining bed. Vertical 
flow is assumed only for the confining bed. Storage 
in the confining bed is considered. Constant dis- 
charge from an aquifer. 


Neuman and Witherspoon (1972) - Design and analy- 
sis of aquifer tests for leaky multiple aquifer sys- 
tems. Observation wells in =quifer and confining bed 
at same distance from pumping well. 


Constant drawdown. 


Jacob and Lohman (1952) - Method uses measure- 
ments of the decreasing flow rate after the well is 


opened. 


Hantush (1959b) - Solutions for drawdown away 
from and discharge at the pumping well for exten- 
sive, circular aquifers. 


Rushton and Rathod (1980) - Analysis by numerical 
methods. 


Variable discharge. 


Hantush (1964b) - Solutions for drawdown corre- 
sponding to three general types of decreasing dis- 
charge. 


Moench (1971) - Convolution integral applied to a 
general discharge function and the method of 
Hantush and Jacob (1955). Representation of dis- 
charge as a step curve and evaluation of the integral 
by summation. 
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Lai and Su (1974) - Includes the effects of storage 
within the pumping well. Presents a general method 
as a convolution integral. Presents methods for 
exponential, pulse function (pumping followed by 
recovery), and periodic (repeated pulse) discharge as 
unevaluated integrals of complex functions. No table 
of values. Presents 19 profiles of the cone of depres- 
sion and 16 curves of drawdown in the pumping well 
as examples of the three discharge functions. 


Nonleaky confined aquifer with homogeneous porous media 
and vertical flow components. 


The pumping well, by the manner of its construction, is 
connected to only a part of the vertical extent of the aquifer. 
Consequently, vertical flow occurs in the vicinity of the 
pumping well. 


Mansur and Dietrich (1965) - Discussion of a series of 
aquifer tests where a fully penetrating pumping well was 
backfilled to create successively smaller partial penetra- 
tions. Analysis of radial-vertical anisotropy through steady 
mined both by electrolytic-analog model and using method 
of Muskat (1946). 


Hantush (1964a) - Application to observation wells piezom- 
eters of the method derived by Hantush (1957). Tables of 
function values. 


Hantush (1961d) - Methods of applying method of Hantush 
(1961a) to analysis of aquifer tests. 


1. Entrance losses. 


Jacob (1947) - Well loss is proportional to the square of the 
discharge. 

Rorabaugh (1953) - Well loss is proportional to the n-th 
power of the discharge. 


Lennox (1966) - Expresses formation loss as a function of 
time through the log approximation to the well function 
(Theis, 1935). 
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2. Inertial effects. 


Cooper and others (1965) - Response to seismic waves. 


Bredehoeft and others (1966) - Inertial and storage effects. 
Ramey (1979), and Kipp (1985). 


Krauss (1974), Shinohara and Ramey (1979), and Kipp 
(1985) - Provides solutions to the oscillatory fluctuations in 
a well after sudden injection or removal of a volume of 
water. 


3. Storage effects. 


Papadopulos and Cooper (1967) - Drawdown in a large- 
diameter pumping well. Storage in the pumping well is an 
important factor in early response. 


II. Unconfined Aquifer 


A. 


Isotropic and homogeneous porous unconfined aquifer 


Boulton (1954a) - A radial, vertical, and time solution as- 
suming all storage at the water table. Most of the discus- 
sion and the limited number of function values are for 
drawdown at the water table. 


Boulton (1954b) - A radial and time method with storage 
throughout the aquifer. A source term in the differential 
equation is referred to as delayed yield or delayed drainage 
at the water table. The delayed yield drainage is the prod- 
uct of an empirical factor with the water-table storage and 
a convolution integral of rate of drawdown and an exponen- 
tial function. 


Boulton (1963) - Same conditions as in Boulton (1954b). 
Type curves and discussion of their use. 


Boulton (1964) - Reply in discussion of Boulton (1963). 
A short table of function values is included. 


Prickett (1965) - Use of the solution of Boulton (1963). Type 
curves and examples. 
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Stallman (1965) - Type curves constructed using analog 
models for two cases where the pumping well is screened 
throughout all of the bottom portion of the aquifer. Storage 
only at the water table and vertical flow. Effects of radial- 
vertical anisotropy are factored into the curves. 


Norris and Fidler (1966) - An example of the use of the 
method of Stallman (1965). 


Anisotropic unconfined aquifer. 


Dagan (1967) - A partially penetrating pumping well. Stor- 
age only at the water table. Anisotropy in the radial- 
vertical plane by a change in scale. 


Neuman (1972) - Fully penetrating pumping well. Storage 
within the aquifer and at the water table. Anisotropy in 
the redial-vertical plane. 


Streltsova (1972) - An interpretation of the a of Boulton 
(1955) as hydraulic conductivity (vertical direction) divided 
by specific yield and a vertical length. 


Neuman (1974) - Partially penetrating pumping well. Stor- 
age within the aquifer and at the water table. Anisotropy 
in the radial-vertical plane. 


Neuman (1975) - Application of the methods of Neuman 
(1972 and 1974). Tables and curves. Interpretation of the a 
of Boulton (1954b and 1963) as not a constant, but varying 
with radial distance from the center of the pumping well. 


Water table in confining layer overlying confined aquifer. 


Cooley (1972) - Interpretation of the a of Boulton (1954b 
and 1963) in terms of properties of an overlying confining 
layer. Numerical models for this situation agree with the 
method of Boulton (1954b). 


Cooley and Case (1973) - The convolution integral of 
Boulton (1954b and 1963) interpreted as the vertical 
velocity at the base of a confining bed with negligible com- 
pressibility. Numerical models indicated that the unsatur- 
ated zone has little effect on flow in the aquifer. 
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Boulton and Streltsova (1975) - Partially penetrating 
pumping well. Storage within the aquifer and at the water 
table but not within the confining bed. Anisotropy in the 
radial-vertical plane. No curves or function values. 
Approximate solutions using method of Boulton (1963). 





III. Other Conditions 


Hantush (1962a) - Flow to a well in a nonleaky confined 
aquifer with a thickness that is an exponential function. 


Hantush and Papadopulos (1962) - How to collect wells 
with lateral (horizontal) screens. 


Bixel and others (1963) - Linear (half-plane) discontinuities 
in hydraulic conductivity or storage coefficient or both. 


Brikowski (1993) - Estimating ground-water exchange 
between ponds or large-scale conduits embedded in uniform 
regional flow. 


Moench and Prickett (1972) - Solutions for estimating 
movement of ground water from a pond or large scale 
radius conduit. 


Zlottnik (1994) - using Dimensional analysis to determine 
and interpret slug test data in anistropic aquifers. 
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Appendix B 


List of Nomenclature and Symbols 


Symbol Dimension 


A L? 
B L 
B l 
C 1 
Do L 
F(B,q) _..... 
H  L 
Hp iL 
H(u,B) ..... 
JQ evens 
 4we 
LT! 
LT" 
K LT? 


Description 
Cross-sectional area 
(Tb/K')1/2 


Coefficient of head loss 
linearly related to the flow 


Coefficient of head loss due to 
turbulent flow in the well, aquifer, 
and across the well screen 


1.5b\ | r(,Kz/Ky) 
F function of B,a 


Change in head in pumping/ 
slug well 


Initial head rise in pumping/ 
slug well 


H function u,B 


Zero order Bessel function of 
the first kind 


First-order Bessel function of 
the first kind 


Hydraulic conductivity of aquifer 


Hydraulic conductivity of confining 


bed 


Hydraulic conductivity of upper 
confining bed 


Equations 
(2) 


(32, 34, 38, 40) 


(5) 


(4) 


(28, 29) 


(28) 


(24, 28) 


(49, 52) 


(24, 27, 54) 


(24, 27) 
(1, 2, 3, 37) 


(34, 36, 37, 40, 
43, 44, 46, 48) 
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Zero-order modified Bessel function 


Degree of anisotropy, equal to K,/K, 
of the second kind 
L (leakance) function of u,v 


Discharge rate 


Hydraulic conductivity of lower 
Storage coefficient 

Storage coefficient 

Storage coefficient of upper 


confining bed 
confining bed 
Specific yield 


LT" Vertical hydraulic conductivity 
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LT! 


L3LT" 


K, LT" Horizontal hydraulic conductivity 


K" 
Ky 

K, 

Ko(x) 
L{u,v) ..... 
Q 

S 

S 

S 

S 

s" 

Sy 
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Ss L* 
Ss L* 


Specific storage of confining beds 
Specific storage of aquifer 


T L*T' Transmissivity 


W (well) function of u 


W (well) function of u for pumped 
control well 


W (well) function of u for image 
control well 


W (well) function of u,r/B 


Zero-order Bessel function of the 
second kind 


First-order Bessel function of the 
second kind 


Aquifer thickness 


Initial saturated thickness of 
unconfined aquifer 


Thickness of confining bed 


Appendix B * * 


-«~ 
* 
pe DS 
oo 
© 
_ 
= 


SPESRSSSES 


—_ pa 
a 
—_ 
> 


BSSSESERE 
TSSEBSERY 


SSSRSERSE 


oo 
=) 
= 
& 


(68) 


(68) 


(32, 38) 


(24, 27) 


(24, 27) 


(1, 4, 37) 


(62, 63, 35, 66) 


(34, 36, 40, 43, 
44, 46, 48) 


(51) 
(51) 
(2, 3) 
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h L 
n L 
r L 
Ti L 
T9 L 
Ic L 
Tj L 
Ip L 
Tw L 
Tw L 
ag L 
8 L 
8 L 
8} L 


Distance, radial, 1 

Distance, radial, 2 

Radius of pumping/slug well casing 
or open hole in the interval where 
water level changes 

Radial distance from center of 
image pumping well 

Radial distance from center of a 
pumping/slug well 


Radius of pumping/slug well screen 
or open hole 


Effective radius of a pumping well 
Drawdown of head 


con 
oo 
~~ 
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A 
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(25, 26, 30, 31) 


(70) 


(69) 


(24, 25, 31) 
(71) 


(6, 9, 11, 13, 14, 
15, 20, 32, 38, 
41, 45, 47, 49, 
54, 59, 65, 66) 


Corrected drawdown, equal to s-(s2/2b) (65) 


Residual drawdown 


(22, 23) 


Dimensionless drawdown equal to 4xTs/Q (59) 





Drawdown in observation well (67, 68) 


é 
s. 


8p L Drawdown component caused by 
pumping control well (67) 


Bw L Drawdown in the pumping/slug well (5) 


t T Time (7, 10, 12, 13, 
14, 16, 20, 21, 
26, 30, 33, 39, 
42, 44, 50, 54, 
60, 61, 69, 70, 
71) 

t T Time since pumping started (23) 

t T Time since pumping ceased (23) 

. Dimensionless time with respect toS (55, 56, 61) 

ty sees Dimensionless time with respect to Sy (60) 

ti T Time, elapsed, 1 (16) 

te T Time elapsed, 2 (16) 

a... —_ r2S/4Tt (6, 7, 8, 10, 12, 
22, 32, 33, 35, 
38, 39, 41, 42, 
44, 49, 50, 52) 

le r2S/4Tt (22) 

-. ‘a Variable of integration (24, 27, 29) 

ugly) —.... Defined in equation (54) (54, 55) 

Unly) —.... Defined in equation (55) (54, 56) 

Up eens r2S/4Tt for pumped control well (69) 

eee r2S/4Tt for image control well (70) 

te r/w(K/Tb )12 (35, 36, 41, 43, 


44) 





pont 
& 


Vv L/T Filux-specified discharge (2) 
fl,v) L/T Average linear velocity (3) 
aa Independent variable in definition 
of erfc(x) (53) 
> ae r(K/Tb )” (45, 46, 47, 48) 
> ee Variable of integration (6, 8, 35, 52, 53, 
54, 55, 56, 57, 
58) 
i: Variable of integration (32) 
. : ae tw Sir," (24, 25, 27, 28, 
29, 31) 
>. : Tt/r,” (24, 26, 28, 29, 
— 30) 
> -~ (r/4bX KS, )” (37) 
a Kpr*/b’ (54, 55, 56) 
> - (r/4X(K S/o TS) + 
(KS / TS)”) (49, 51, 52) 
a r’K,/(K;b’) (63) 
 - Root of equation (56) (55, 57) 
~~ Root of equation (57) (56, 58) 
a ee Change in parameter (14, 15, 18, 19, 
(finite difference) 23) 
Ast L Change in drawdown over one log 
cycle of time (18) 
Asr L Change in drawdown over one log 
cycle of radial distance (19) 
Au)... Function of u defined in equation (26),(24, 27, 29) 
aa S/Sy (56, 56, 57, 58) 





Appendix C 
Glossary of Geohydrologic Concepts 
and Terms 


The following definitions and concepts of geohydrologic terms are princi- 
pally from Fetter (1994). The book by Fetter, Applied Hydrogeology 
(1994), is also an excellent source for ground water analysis and inter- 


pretation. 


Anisotropy - Anisotropy is that condition in which significant proper- 
ties are a function of direction. Anisotropy is common in sedimentary 
sequences in which hydraulic conductivity perpendicular to the bedding 
planes is less than the hydraulic conductivity parallel to the bedding. 


Aquifer - An aquifer is a saturated geologic unit that has sufficient 
permeability to transmit water at a substantial rate. An aquifer is com- 
monly defined, in terms of water yielding capacity, as a formation, group 
of formations, or part of a formation that contains sufficient saturated 
material to yield significant quantities of water to a well or 


springs 


Aquifer, leaky - A misnomer, but used here and in aquifer test litera- 
ture to refer to a confined aquifer that receives leakage from adjacent 
confining beds when the aquifer is stressed by a pumping well. (See 
confining bed, leaky.) 

Aquitard - See preferred term, confining bed or leaky confining bed. 
Artesian - Artesian is synonymous with confined; artesian aquifer is 
equivalent to confined aquifer. An artesian well is a well deriving its 
water from an artesian or confined aquifer. The water level in an arte- 
sian well stands above the top of the artesian or confined aquifer it 
penetrates. 


Confining bed - A confining bed is a geologic unit with minimal perme- 
ability. These beds are not permeable enough to yield significant quanti- 
ties of water to wells or springs. The permeability of aquifers and con- 
fining beds is not precisely defined in a quantitative sense, but a confin- 
ing bed has distinctly less permeability than the aquifer it confines. 
Other terms that have been used for beds with minimal permeability 
include aquitard, aquifuge, and aquiclude. 
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Confining bed, leaky - A leaky confining bed yields a significant quan- 
tity of water to the adjacent aquifer when the aquifer is stressed by a 
pumping well. 


Drawdown - Drawdown is the difference between the static water level 
and the water level after pumping has begun. 


Effective radius - The effective radius of a well is that distance, mea- 
sured radially from the axis of the well, at which the theoretical draw- 
down based on the logarithmic head distribution equals the actual draw- 
down just outside the well screen (Jacob, 1947). From the time intercept 
of the time drawdown logarithmic plot with the zero drawdown line, the 
effective radius of the pumping well can be determined by the following 
equation from Jacob (1947, equation 25, p. 1059): 


ry? = 2.25\f{Tt,S) 
Equilibrium, state of - See flow, steady. 


Flow, steady - Steady flow occurs when, at any point, the magnitude 
and direction of the specific discharge are constant in time. 


Flow, unsteady - Unsteady or nonsteady flow occurs when, at any 
point, the magnitude or direction of the specific discharge changes with 
time. 


Ground water, confined - Confined or artesian ground water is under 
pressure significantly greater than atmospheric, and its upper boundary 
is the bottom of a bed of distinctly lower hydraulic conductivity than 
that of the bed in which the confined water occurs. 


Head, total - The total head of a liquid at a given point is the sum of 
three components: (1) elevation head, which is equal to the elevation of 
the point above a datum; (2) pressure head, which is the height of a 
column of water that can be supported by the static pressure at the 
point; and (3) the velocity head, which is the height the kinetic energy of 
the liquid is capable of lifting the liquid. 


Homogeneity - Homogeneity is synonymous with uniformity. A mate- 
rial is homogeneous if its hydrologic properties are identical everywhere. 
Although no known aquifer or confining bed is homogeneous in detail, 
models based on the assumption of homogeneity have been determined 
empirically to be valuable tools for predicting the approximate relation 
between ground-water flow and hydraulic head in many flow systems. 
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Hydraulic conductivity - The hydraulic conductivity of a medium is 
the volume of water at the existing kinematic viscosity and density that 
will move in unit time under unit hydraulic gradient through a unit area 
measured at right angles to the direction of flow. Hydraulic conductivity 
has dimensions of velocity. 


Isotropy - Isotropy is that condition in which all significant properties 
are independent of direction. 


Nonequilibrium, state of - See flow, unsteady. 


Observation well - An observation well is open to the aquifer through- 
out a given vertical distance. The water level in an observation well 
reflects the average head in the aquifer profile that is occupied by screen 
or perforated casing (Hantush, 1961a). 


Permeability, intrinsic - Intrinsic permeability is a measure of the 
relative ease with which a porous medium can transmit a fluid under a 
potential gradient. It is a property of the medium alone and is theoreti- 
cally independent of the nature of the fluid and of the force field causing 
movement. 


Piezometer - A piezometer is a small-diameter pipe open to the aquifer 
only at its lower end (Hantush, 1961a). 


Porosity, effective - Effective porosity is the amount of interconnected 
pore space available for fluid transmission. 


Potentiometric surface - The potentiometric surface at a point is 
defired by the level to which water will rise in a tightly cased well or 
piezometer. 


Pumping/slug well - The pumping/slug well of an aquifer test is the 
well through which the aquifer is stressed, for example, by pumping, 
injection, or change of head. 


Saturated zone - The saturated zone is a zone beneath the ground 
surface in which all voids, large and small, are ideally filled with water 
under pressure greater than atmospheric. 


Specific capacity - The specific capacity of a well is the discharge per 
unit drawdown. The specific capacity usually decreases both with time 


and discharge. 
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Specific discharge - Specific discharge is the rate of discharge of 
ground water per unit area of porous medium measured at right angles 
to the direction of flow. 


Specific storage - The specific storage of a confined aquifer is the vol- 
ume of water released from or taken into storage per unit volume of the 


porous medium per unit change in head. 


Specific yield - The specific yield of a rock is the ratio of the volume of 
water that the saturated rock will yield by gravity to the volume of the 
rock. Specific yield is determined by tests of unconfined aquifers and is 
the change that occurs in the volume of water in storage per unit area of 
unconfined aquifer as the result of a unit change in head. Such a change 
in storage is produced by the draining or filling of pore space and is, 
therefore, dependent on particle size, rate of change of the water table, 
time, and other variables. 


Storage coefficient - The storage coefficient is the volume of water an 
aquifer releases or takes into storage per unit surface area of the aquifer 
per unit change in head. In a confined aquifer, the water derived from 
storage with decline in head comes from expansion of the water and 
compression of the aquifer. In an unconfined aquifer, the volume of water 
derived from or added to the aquifer by these processes is much smaller 
compared to that involved in gravity drainage. 


Storativity - Synonymous with storage coefficient. 


Transmissivity - Transmissivity is the rate at which water of prevailing 
kinematic viscosity is transmitted through a unit width of the aquifer 
under a unit hydraulic gradient. 


Unsaturated zone - The unsaturated zone is the zone in which water is 
under less than atmospheric pressure. This zone is also referred to as 
the vadose zone and the zone of aeration. 


Vadose zone - See preferred term, unsaturated zone. 


Water table - The water table is that surface in an unconfined aquifer 

at which the water pressure is atmospheric. It is defined by the level at 
which water stands in a well that penetrates the aquifer just far enough 
to hold standing water. 


Well loss - A component of drawdown in a discharging well. Well loss is 
the loss of head in a pumping well due to turbulent flow that accompa- 
nies the flow of water through the aquifer, screen, and upward inside the 


casing to the pump intake. 
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